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ABSTRACT 

POLARIZATION  OF  THE  COSMIC  BLACKBODY  RADIATION 

AT  3.2  CM. 

GEORGE  PETER  NANOS,  JR. 

* 

A  3.2  Dicke  Radiometer  configured  as  a  polarimeter  was 
used  to  make  measurements  of  the  linear  polarization  of  the 
Primeval  Fireball  along  a  declination  of  ^+0.35  N.  The 
resolution  of  the  instrument  was  15  in  declination  and  1  hour 
in  right  ascension.   Harmonic  analysis  of  the  Stokes  parameters 
of  the  radiation,  Q  and  U,  sets  limits  on  sidereal  variations, 
which  in  turn  limit  the  possible  contributions  of  cosmological 
origin.   Results  are  related  to  the  isotropy  of  the  cosmic 
background  radiation  and  the  symmetry  of  the  universe. 


I.   INTRODUCTION 

1.1  HISTOKY  OF  THE  FIREBALL 

In  1929  it  became  clear  from  the  work  of  Edwin  Hubble  that  the 
universe  is  expanding,  and  that  the  recession  velocity  of  a  galaxy 
is  nearly  a  linear  function  of  its  distance.   This  result  gave 
impetus  to  the  development  of  expanding  universe  cosmologies  based 
on  solutions  to  Einstein's  equations.   Of  the  many  ensuing  models, 
among  the  most  successful  and  well  established  is  the  homogeneous 

and  isotropic  universe  with  cosmological  constant  equal  to  zero, 

2  3 

based  on  the  work  of  Friedman  and  Lemaitre.  ' 

One  important  feature  of  this  model  is  the  existence  of  a 
singularity  at  t=0;  i.e.,  at  some  finite  time  in  the  past  the  scale 
factor  in  the  universe  was  small  and  the  density  of  matter  and 
radiation  was  very  high.  This  feature  is  known  as  the  "fireball" 
hypothesis  and  the  subsequent  expansion  and  evolution  of  the 
universe  is  called  the  "Big  Bang". 

During  the  early  stages  of  the  fireball  expansion,  the  matter 
in  the  universe  was  extremely  hot,  much  greater  than  10    K,  and 
was  in  thermal  equilibrium  with  the  radiation.   Being  thermal,  this 
radiation  had  a  Planckian  spectrum  given  by 

8nhv13 

dE.    =  3 hTkT-  dV  (1-D 

1  c  ell 

The  first  description  of  this  type  of  blackbody-radiation  filled 

k 
model  was  given  by  Tolman  in  1931.   Later  he  showed  that  under 

adiabatic  expansion  the  radiation  spectrum  is  still  given  by  (l.l), 


where  the  temperature  T  is  now  given  by 

a         . 

T'(t)  =  T  -TTT  =  T  —-   ,  (1.2) 

x  '  o  a(t)    o  1+z  '  x 

where  a(t)  is  the  scale  factor  of  the  universe  at  time  t,  and  1+z 

is  the  redshift  due  to  the  first  order  Doppler  effect  from  the 

expansion,  viewed  by  a  cornoving  observer. 

Alpher  and  Gamow  in  19^8  needed  the  fireball  in  their  attempt 

to  explain  element  formation  in  the  early  universe  and  in  their  theory 

derived  a  temperature  of  5  K  for  the  remnants  observed  today,  a  figure 

o  5 

quite  close  to  the  2.7  K  obtained  experimentally.    Interest  in  this 

result  waned  as  other  explanations  of  element  formation  proved  more 

attractive  until,  working  independently,  R.  H.  Dicke  and  his  colleagues 

at  Princeton  were  drawn  to  the  idea  to  explain  the  reprocessing  of 

heavy  elements  into  hydrogen  in  a  closed  oscillating  universe. 

At  the  same  time  that  two  of  these  co-workers,  P.  G.  Roll  and 

D.  T.  Wilkinson,  were  building  a  radiometer  to  search  for  this 

radiation,  Penzias  and  Wilson  at  Bell  Laboratories  realized  that 

what  they  previously  had  considered  to  be  excess  receiver  noise  at 

7.35  cm.  and  had  tried  to  eliminate,  was  in  fact  the  remnant  of  the 

7 
primeval  fireball  that  the  Princeton  group  was  searching  for. 

Soon  afterward  the  3.2  cm.  measurement  of  Roll  and  Wilkinson  con- 
firmed the  result  and  showed  that  the  spectrum  over  this  range  was 

consistent  with  that  of  a  2.7  K  blackbody.   Following  this,  a  host 

9-28 
of  other  measurements  have  been  performed/    all  of  v/hich  tend  to 

support  the  hypothesis  to  a  greater  or  lesser  degree,  at  least  in 

the  Rayleigh- Jeans  part  of  the  spectrum.  At  wavelengths  greater  than 


a  few  tens  of  centimeters  or  shorter  than  a  millimeter,  the  situation  is 
not  as  clear.  At  longer  wavelengths  contributions  due  to  sources 
within  the  galaxy  make  separating  out  the  blackbody  difficult  while 
short  wavelength  observations  suffer  greatly  from  absorption  and 
emission  in  the  atmosphere.   If  one  excludes  the  rocket  data  at  very 
short  wavelengths,  the  reliability  of  which  seems  open  to  question 

at  the  present  time,  Peebles  has  shown  that  the  available  data  supports 

29 
the  fireball  reasonably  well. 

ANISOTROPY 

Accepting  the  idea  that  one  is  looking  at  fossil  radiation  from 
early  in  the  history  of  the  universe,  it  is  natural  to  ask  what 
information  it  carried  from  this  early  epoch.  As  the  early  universe 
expanded  and  cooled,  at  one  point  the  primeval  plasma  recombined  and 
the  mean  free  path  for  photon  scattering  became  very  large,  so 
large  in  fact  that  most  of  the  fireball  photons  we  observe  were 
probably  last  scattered  during  the  recombination  at  z=1000.   Only 
a  subsequent  reionization  of  the  matter  could  bring  the  scattering 
epoch  closer.   Therefore,  the  spatial  distribution  reflects  the 
distribution  and  lemperature  of  matter  at  this  early  time. 

Tests  for  isotropy  in  the  distribution  of  the  radiation  have 
been  conducted  by  several  investigators  including  the  Princeton 
group.      At  the  present  time  no  significant  departure  from 
spherical  symmetry  has  been  found  for  angular  resolutions  on  the 
order  of  1  of  arc.  In  addition  to  early  symmetry  of  the  fireball, 
the  data  can  be  used  to  deduce  our  peculiar  velocity  relative  to  the 

local  cornoving  frame.  This  effect  appears  as  a  2k   hr.  asymmetry  with 

39-Ul 
a  hot  spot  in  the  direction  of  motion  due  to  the  added  Doppler  shift. 


1.2  COSMOLOGICAL  CAUSES  OF  POLARIZATION 

As  was  pointed  out  by  Martin  Rees  in  i960,  "  the  type  of  expansion 
asymmetries  which  would  give  rise  to  anisotropics  in  the  microwave  back- 
ground would  also  give  rise  to  linear  polarization  of  the  radiation 
through  the  mechanism  of  Thomson  scattering.  Following  his  lead 
and  using  much  the  same  notation,  I  will  briefly  demonstrate  the 
existence  of  a  linearly  polarized  component  pf  the  background,  given 
an  asymmetric  expansion. 

In  this  I  will  consider  only  one  particular  class  of  asymmetric 

universes,  those  with  axially  symmetric  Euclidean  metrics  of  the 

2     2       2    2     2        2   2 
form  ds  =  dt   -  A(t)   (dx  +  dy  )  -  W(t)   dz  and  containing  no 

magnetic  field.   Models  of  this  type  have  been  investigated  extensively 
by  Thorne  (1967),   Jacobs  (1968),   and  others  and  are  particularly 
attractive  because  of  their  simplicity,  their  close  relation  to  the 
standard  Einstein-de  Sitter  universe,  and  the  straightforward  way 
in  which  they  lead  to  a  polarization.   In  the  case  of  three 
inequivalent  directions,  the  situation  would  be  much  more  complex 
(there  being  azimuthai  dependence  to  the  equations)  and  the  calcula- 
tions correspondingly  more  difficult. 

Taking  our  given  metric  we  can  first  define  Hubble ' s  constants  in 
each  direction,  a  =  dA/Adt  and  w  =  dw/v/dt.  From  these  we  can  compute 
an  average  expansion  constant  h  =  (2a+w)/3  and  a  differential  constant 
Ah  =  (w-a)  and  using  them,  can  discuss  the  existence  of  a  temperature 
asymmetry  in  this  type  of  universe.  As  v;as  mentioned  in  the  preceding 
section,  Tolman  showed  that  in  the  Einstein-de  Sitter  case  the 
temperature  of  radiation  is  inversely  proportional  to  the  scale 


factor,  which  in  turn  is  proportional  to  v  '  .   In  our  case  V  the 

2 
volume  is  proportional  to  A  W  and  in  a  similar  way  we  can  write 

T  a  (A^)"  ""A  (1.3) 

where  T  is  an  average  temperature  over  all  directions  of  observation. 
At  any  time  t,  one  can  define  an  approximate  isotropy  parameter  equal 
to  AT/T  =  e.   AT  is  the  temperature  difference  between  the  A  and  W 
directions  and  can  be  written  ! 


AT  =  t  d/dt  (T  -  T  ) 
B  '    v  w   a.' 

a   t  d/dt  (l/W  -  l/A)  (1'4) 

s 

~  t  Ah/tA^)1/3   . 
—   s 

Here  t  is  the  time  since  the  last  scattering,  and  is  given  by  \    , 

the  mean  free  path  of  the  photons,  divided  by  the  velocity  of  light; 

i.e. 

i 

t 
s    c 

1 


(1.5) 


nl(t)aTc 

where    n  =  the  number  density  of  electrons 

l(t)  =  the  percentage  of  them  which  are  ionized 
ff_  =  the  Thomson  scattering  cross  section. 

Putting  this  form  into  the  previous  equation,  we  have  finally, 

e  =  (aS/)1/3  A  T 

Ah  ^-6) 

-  ni(t)oTc   * 

To  obtain  a  more  precise  result  requires  consideration  of  the 
cosmological  redshift  in  an  asymmetrically  expanding  universe  as  is 
done  in  Appendix  A2.  In  particular  for  the  axial  model  being 
considered,  the  temperature  of  the  radiation  as  a  function  of  the 


direction  of  observation  'with  respect  to  the  axis  of  the  universe  is 

shown  to  be 

2  2 

A  W       ^i 

T  =  T  [-f-  sin2  6  +  -§-  cos2e]"  2  ,  (1.7) 

A  W 

in  the  absence  of  scattering. 

Throughout  this  discussion  the  zero  subscript  will  refer  to 

the  value  of  a  quantity  at  the  present  time.   In  the  event  that  the 

expansion  asymmetry  is  small,  as  it  is  in  the  case  we  are  considering, 

this  equation  can  be  cast  in  a  little  more  convenient  form  as 

T  =  T   (1  +  e  sin2  9).  (1.8) 

Taking  this  type  of  temperature  distribution  as  a  starting  point, 

the  polarization  of  the  microwave  background  through  the  mechanism  of 

Thomson  scattering  can  now  be  demonstrated.  As  shown  in  Figure  1.1, 

we  will  define  the  polarization  of  an  electromagnetic  wave  from  the 

microwave  background  as  seen  by  an  observer  in  terms  of  the  plane 

containing  the  line  of  sight  to  the  observer  and  the  direction  of 

asymmetric  expansion.  The  direction  perpendicular  to  this  plane  is 

called  the  'a'  direction  while  the  direction  which  lies  in  this  plane 

and  is, of  course,  perpendicular  to  the  direction  of  propagation  is 

called  the  coplanar  or  'w'  direction.   Following  P.ees,  I  will 

write  the  temperature  of  the  radiation  polarized  in  each  direction  in 

the  following  way: 

T  =  T  (1  +  e  sin2  9) 
a     v     a       ' 

Tw  =  T  (1  +  ew  sin2  9)  ,  (1.9) 

where  T  now  refers  to  T 


w 


A 


-\    -    A 


FIG. 4.1 


DEFINITION  OF  THE  POLARIZATION  REFERENCE 
PLANE 


By  doing  so, I  can  characterize  the  polarization  completely  by  means  of 
the  parameters   e   and  e   .   In  fact,  the  evolution  of  the  polarization 

3.  W 

with  time  may  be  ascertained  by  constructing  and  solving  the  differen- 
tial equation,  describing  the  change  of  these  parameters  with  expansion 
and  scattering  of  radiation.   It  must  be  understood  that  this  type  of 
approach  is  only  valid  when  the  temperature  is  low  enough,  say 

Q  o 

10   K,  for  elastic  Thomson  scattering  to  be  the  primary  mode  of 
interaction  with  the  photons  and  when  the  asymmetry  in  the  expansion 
is  small  enough  that  the  above  expansions  of  the  temperature  distribu- 
tion are  valid.   In  light  of  the  small  observed  temperature  asymmetry 
this  latter  condition  doesn't  seem  too  unreasonable.  At  this  time  there 

doesn't  seem  to  be  any  evidence  for  a  very  hot  plasma  greater  than 

9  o 
T  =  10   K  at  a  recent  epoch. 

If  we  express  the  temperature  of  the  radiation  in  the  two 

polarizations  of  the  microwave  background  in  the  way  described  above, 

we  can  compute  the  effect  on  the  asymmetry  parameters  e  and  e   of 

a  single  Thompson  scattering.  In  Appendix  A3  the  result  in  matrix 

form  after  integration  over  all  angles  of  incidence  -is  shown  to  be 

0   0 


(:■)  ■  C 


(1.10) 


-1/2  7/30/    \e 

where  the  primed  quantities  refer  to  the  values  after  scattering. 

The  most  notable  effect  in  this  calculation  is  that  the  asymmetry 

in  the  orthogonal  'a'  component  is  destroyed  by  the  scattering  process 

while  the  asymmetry  parameter  in  the  coplanar  polarization  is  preserved. 

If  to  begin  with,  both  parameters  had  been  equal  to  e  ,  the  effect 

of  the  scattering  would  have  been  to  increase  the  maximum  polarization 


from  zero  to  l/30  e,  which  demonstrates  the  effect  of  scattering  in 
producing  a  polarization. 

Because  we  know  the  time  scale  over  which  this  change  occurs, 
we  can  write  an  expression  for  the  rate  of  change  of  the  asymmetry 
parameters  with  time  as  f ollows : 


ft). 


d_ 

dt  \~w/     At 


ft) 


t 
s 


'e  ' 
a 


e  • 

»  w 


ft) 


=  nl(t)  aTc  (-1/2  -23/30]  (^ J.  (1.11) 


In  addition  to  the  effect  of  scattering  on  e   and  e  ,  one  must 

a      w ' 

consider  the  effect  of  expansion  even  though  this  is  not  what  leads 
to  the  existence  of  a  polarization.   From  equations  1.5  and  1.6  we 
find  that  the  growth  of  a  temperature  asymmetry  with  time  due  solely 
to  expansion  is  t  Ah.  This  implies 


dt   IsJexp  ■  -4r^ 


At 
1 


m     tsAh  (j)   m     Ah  (i}  >  (1>12) 

t 
s 

Combining  this  with  the  previous  result  we  obtain  a  final  differential 

equation  for  e   and  e  : 
a      w 

h     ft)  ■  ffl  4)  *  "-^^   (-1/2  -23/30)  (<,)•  »■«> 

Using  the  result  from  Appendix  Al  that  All  =Ah./A  W  and  noting  that 
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n  =  n_/A  w  simplifies  the  expression  to 


0 

d_ 
dt 


•W     ^o    1     V(t)CTTC     A1    °     \  /ea\ 

(<)■■£  a*-^-  (-1/2-23/30 j  («: )  •  <-*> 


It  is  also  shown  in  Appendix  Al  that  A  and  W  can  be  expressed  in 
the  following  way: 

A  =  A   (t/t  )2/3 

00  i  (1.15) 

-1/2 
W  =  A  +  qA  I     , 

which  demonstrates  their  time  dependence  explicitly,  leaving  only  one 

unknown  function  of  time  in  the  equation  l(t).  The  use  of  q  here 

should  not  be  confused  with  the  deceleration  parameter.   I  only  use 

it  for  consistency  with  Thome's  notation.  Since  we  have  no 

theoretical  or  experimental  handle  on  the  functional  dependence  of 

l(z),  I  will  follow  Rees  and  set  it  equal  to  one  in  order  to  obtain  a 

useful  approximate  solution.  Adding  this  to  the  above,  we  obtain  the 

following  two  differential  equations  for  the  asymmetry  parameters: 


de  n  arnc  Ah 

a       o  I o 

dt  A3    (t/t    )2+qA  3/2(t/t    )      ^      "  A3    (t/t    )2+qA   3/2(t/t    ) 

o     '    o  o  '    o  o     '    o  O  '    o 


dew  noQTc ,  1  23       ,       ^o 


dt      U  A3    (t/t    )2+qA   3/2(t/t    )      '    2      ?a  +   30   £^     '   A3    (t/t    )2+qA   3/2(t/t 
0N   '    o  o        v   '    o  o     '    o  o  ' 

(1.16) 


Solutions  to  these  are  easy  to  obtain  and  are  given  by 
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All  /  at  \      n  amct   /a A  3'2 


e     =  — 1  + 

noV  \  A3/2t 

o  ' 


3/2 


15Ah  /  at  \     n  act    ,  A 

0  I    -,  °  \       o  T      o/q  o  /-,    nr,\ 

*.  =  ! 1  +    .  3/2  t     1  '  ^-17) 


*oV  V  Ao 


which  asymptotically  approach 
Ah 


e 
a 


o  T 


(1.18) 


15Ah 
e  t  =  o 

7nooTc 

o  /p 
In  the  event  that  q_  is  a  small  number,  as  long  as  t  >  qt  /A  '  , 

these  values  are  reached  and,  in  fact,  it  is  shown  in  Appendix  A2 

that   e  being  small  is  equivalent  to  requiring  that  q  be  small. 

From  the  asymptotic  solutions  we  can  now  compute  e  =  -g"  (e  +  e  ) 

and  the  degree  of  polarization  (T  -  T  )/l. 

-1-  v  a   w' ' 

11    o 
e 


7  noaTc 


(1.19) 
o  Ah 

(T  -  T  )/T  =    S sin  0 

v  w   a"       in  amc 

o  T 

8  I  I   •  2  „ 

=   II  lel  sin  e 

This  is  the  scattering  induced  polarization  which  we  wished  to 
demonstrate.  Although  the  assumption  that  l(t)  is  identically  equal  to 
one  is  unrealistic,  the  solution  so  obtained  may  now  be  applied  to  more 
tenable  pictures  of  the  evolution  of  the  universe  to  test  for  an 
anisotropic  expansion. 
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II .   APPARATUS 

2.1  THE  MEASUREMENT  OF  LINEARLY  POLARIZED  RADIATION 

In  order  to  measure  and  understand  the  effect  of  a  cosmological 
linear  polarization,  we  must  first  understand  how  polarized  radiation 
is  measured.   The  simplest  kind  of  polarimeter  is  shown  in  Figure  2.1a. 
It  consists  of  a  polarization  sensitive  antenna,  of  which  the  simplest 
is  the  center  fed,  linear  dipole  shown,  followed  by  a  radio  receiver. 
To  measure  polarization,  you  aim  the  polarimeter  at  a  source  and  rotate 
it  until  you  obtain  the  maximum  signal.   The  angle  between  the  axis  of 
the  dipole  and  some  suitable  reference  position  is  the  angle  of  the 
polarized  signal;  and  the  difference  between  this  signal  and  lowest 
signal  received,  as  the  device  is  rotated  through  360  degrees,  is 
proportional  to  the  magnitude  of  the  polarization.   This  method  is 
fine  and, because  of  its  simplicity, is  probably  the  best  as  long  as 
the  signal  is  strong  enough  to  give  a  clear  difference  between  maximum 
and  minimum. 

In  many  practical  applications  the  signal  is  quite  small  and  has 
to  compete  with  a  large  amount  of  system  noise  from  the  receiving 
equipment  itself.   In  this  case  the  apparatus  in  Figure  2.1b  is 
better.  It  consists  of  two  antennas  of  the  dipole  type  with  their 
polarization  sensitive  directions  oriented  at  right  angles  to  each 
other  and  each  connected  to  its  own  receiver.  When  this  is  aimed  at 
a  source,  a  signal  S.  will  appear  at  the  output  of  each  receiver: 

S1=c(Tu+Tpcos2e+TRE      ) 

(2.1) 

S0  =   c(T     +  T     sin29  h    Trt?_    ), 
2  v   u         p  REC?y' 
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where  T  is  the  temperature  of  the  unpolarized  radiation  from  the  source, 
T  is  the  temperature  of  the  polarized  component, 
G  is  the  angle  of  the  polarization  with  respect  to  antenna  one, 
C  is  a  proportionality  constant  depending  on  the  gain  of  the 

receivers,  the  method  of  detection, ' and  the  receiver  band- 
widths  . 
T^^,  is  the  equivalent  noise  temperature  of  each  receiver. 
Subtracting  the  two  signals  gives  an  expression  for  T  and  9, 

S±-   S£  =  c[Tp  COS26  +  (T^-  TREC2)].  (2.2.) 

If  the  noise  figures  of  the  two  receivers  are  the  same,  the  second 

expression  in  the  brackets  has  mean  zero  and,  for  sampling  times  long 

compared  to  the  reciprocal  of  the  receiver  bandwidth,  a  Gaussian 

distribution.   Thus  the  time  average  of  this  term  is  zero  and  we  can 

define 

Q  =  1     °2  Time  =  T  cos29.  (2.3) 

c  p 

Rotating  the  antenna  system  by  forty-five  degrees  gives  another 

pair  of  signals, 

S  '  =  c[T  +  T  cos2(6  +  rt/k)   +   !__  ] 

1  up  v  '     '  REC 

1  (2.4) 

S2'   =   c[Tu+   Tp   sin2(6  +  Tt/if)   +  TREC    ], 

from  which  a  second  parameter  can  be  defined  in  a  similar  way 

U  =  <  S    '-  S    '   >     =  T     sin29.  (2.5) 

1  2  p  \      si 

c 
In  terms  of  these  parameters  the  polarization  temperature  and  angle 
are  given  as 
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TWO  TYPES   OF  POLARIMETER 
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T   =  /~2        2 
p    Q  +  U 

1        U  <2-6) 

9   =  I   ARCTAN  (^)  . 

It  is  clear  that  these  two  parameters,  known  as  the  Stokes 

parameters,  uniquely  determine  the  polarization -of  a  source  and  thus 

together  form  a  complete  set  of  orthogonal  measurements  to  determine 

this  quantity. 

i 

For  many  applications  this  type  of  polarimeter  works  quite  well, 

but  it  too  has  limitations.   Even  when  one  finds  two  receivers  with 

well  matched  noise  figures,  drifts  in  the  gain  of  the  devices,  which 

are  not  necessarily  random  in  nature,  can  cause  problems.   If  the 

signal  one  is  seeking  to  measure  is  of  the  same  order  as  the  product  of 

the  receiver  noise  and  the  size  of  the  average  fluctuations  in  gain 

between  the  two  receivers,  serious  errors  in  the  measurement  result. 

These  errors  cannot  be  averaged  out.  A  solution  to  this  problem, 

diagrammed  in  Figure  2.2,  is  the  replacement  of  the  two  receivers 

with  a  single  Dicke  receiver  which  is  set  to  switch  synchronously 

between  the  two  antennas.    The  power  entering  each  antenna  is 

written  as 

2„ 


T,  =  T  +  T  cos 

1  u    p 

2 

T0  =  T  +  T  sin  i 

2  u    p 


(2.7) 


The  signal  after  the  switch  and  before  the  receiver  is  given  by 

Q  =  T--  Tg  =  T  cos20  f(cut),  (2.8) 

where  f(ujt)  is  a  periodic  function  at  the  switching  frequency,  with 
amplitude  one,  which  is  dependent  upon  the  switching  characteristics  of 
the  receiver  and  the  modulating  signal.  After  passing  through  the 
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receiver,  the  signal  picks  up  a  contribution  due  to  receiver  noise,  which 
is  a  sum  of  contributions  at  all  frequencies  up  to  the  cutoff  of  the  IF 
and  video  amplifiers.   In  particular,  the  contributions  due  to  gain 
fluctuations  are  at  low  frequencies,  which  means  that  by  choosing  the 
switching  frequency  high  enough,  we  can  isolate  the  wanted  signal  from 
the  gain  fluctuations.  How  if  the  signal  is  processed  through  a  phase 
sensitive  detector  like  a  lock-in  amplifier,  only  the  wanted  signal 
and  those  noise  contributions  with  frequencies  in  a  limited  bandwidth 
around  the  switching  frequency  are  recovered  and  the  contributions 
from  gain  fluctuations  are  excluded: 

S  =  c[Q  +  TREC  (Av)].  (2.9) 

Again  taking  a  time  average  we  can  recover  the  parameter  Q, 

<S>      ■  ,     N 

Q  =    Time  .  (2.10) 

c 
If  the  antennas  are  rotated  by  ^5  degrees,  U  can  be  measured  in  a 
similar  way. 

The  above  progress  is  not  gained  without  cost.   By  switching  between 
the  two  antennas  we  are  only  looking  at  each  one  half  of  the  time  and 
effectively  throwing  away  half  the  available  signal.  The  minimum 

detectable  signal  in  a  radio  receiver  as  a  function  of  system  noise 

'(6 
temperature,  integration  time,  and  bandwidth  is  given  by 

T 

AT.   =  REC     .  (2.11) 

mm    /.  .. —  v 

/AvAt 

This  tells  us  that  if  we  have  a  receiver  with  bandwidth  Av  and  we  wish 

to  measure  Q  or  U  to  an  accuracy  AT  then  we  must  carry  out  the 

above  time  average  for  a  period  At  in  order  to  eliminate  the  effects 

of  TR  .   For  the  Dicke  receiver,  since  we  are  throwing  away  half  the 
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signal,  things  became  worse  by  a  factor  of  two: 

T 

AT.   =  2  EEC     .  (2.12) 

111111     n — tz 

/AvAt 

If,  instead  of  a  wide  band  pre -multiplication  amplifier,  a  narrow  band 

tuned  amplifier  is  used,  then  the  first  harmonic  at  the  switching 

frequency  alone  is  extracted  and  the  corresponding  loss  in  sensitivity 

is  k/rr  /2   .   This  loss  is  more  than  made  up  for  by  the  freedom  from 

overloading  and  simpler  amplifiers  made  possible  by  this  approach. 

Finally,  if  the  switching  is  also  done  sinusoidally  (something  which 

is  often  done  as  a  matter  of  convenience),  the  sensitivity  is  further 

reduced  by  the  factor  U/n.  Putting  these  together  gives  the  expression 

for  a  sinusoidally  switched  Dicke  radiometer  with  narrow  tuned  video 

amplifier 

T 
AT.   =  2f2~       REC      .  (2.13) 

mn         /A^At- 
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2.2  TERRESTRIAL  OBSERVATION  OF  A  COSMOLOGICAL  LINEAR  POLARIZATION 

Next  we  ask,  what  kind  of  signal  does  polarized  background 
radiation  produce  in  an  earthbound  radiometer,  given  the  orientation 
of  the  axis  of  the  asymmetric  expansion  with  respect  to  the  rotation 
axis  of  the  earth,  the  latitude  of  the  observer,  and  the  direction 
in  which  the  observer  is  looking.   In  order  to  simplify  the  question 
and  the  subsequent  analysis,  I  will  assume  that  the  earthbound 
observer  will  make  all  his  observations  by  pointing  his  polarimeter 
at  the  zenith.   It  will  turn  out  later  that  there  are  sound  experimental 
reasons  for  this  procedure  in  addition  to  the  computational  ones. 

If  we  were  to  step  away  from  the  earth  for  a  moment  and  aim  our 
polarimeter  in  a  direction  at  an  angle  9  with  respect  to  the  axis 
of  the  universe,  and  if  we  aligned  our  polarimeter  with  the  sensitive 
axis  of  the  antenna  aligned  with  the  coplanar  direction,  then  from 
equations  1.9  we  find  that  we  would  measure  a  Stokes  parameter  Q, 
given  by 

Q  =  T  -  T 

W   a  (2.H0 

p 

=  T(e  -  e  )  sin  9, 
v  w   a       ' 

or  in  an  equivalent,  and  for  my  purposes  more  convenient,  form  as 

T  -  T  =  (T  -  T  )    sin  9.  (2.15) 

w   a   v  w   a'max  v    ' 

Nov;,  stepping  back  on  the  earth  which  is  rotating  about  an  axis,  which 

forms  an  angle  $  with  respect  to  the  symmetry  axis  of  the  universe,  I 

observe  the  zenith  which  is  a  direction  whose  declination  is  given  by 

the  latitude  of  the  point  of  observation  and  whose  right  ascension  is 

just  the  local  sidereal  time  at  the  observation  site.   Figure  2.3  is  a 

schematic  representation  of  the  situation,  showing  the  experimental 

geometry. 
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Given  this  situation,  I  show  in  Appendix  Ak   that  the  polarimeter 

will  measure  a  quantity 

T,-  T„  =  (T  -  T  )     |[cos26(l  -   |  sin2§)  +  sin2§(l  -  \     cos26)cos2(c<  -  t' 
1   2   v  w   a'max  i      x     2  v     2      '  v  o 

2  2 

-  l/2  sin2$  sin25   003(0-   -   t)]   cos    {3   -    [sin   $  sin  6    sin2  (c*   -   t) 

2  / 

-  sin2§  cos  6  sin(a  -  t)]  sin  f3  }  , 

°  >  (2.16) 

where  t  is  the  local  mean  sidereal  time  at  the  observation  site, 

p  is  the  angle  between  the  sensitive  direction  of  the  antenna 

and  the  direction  of  the  North  Celestial  pole, 
6  is  the  declination  of  the  zenith,  and 
a     is  the  right  ascension  of  the  symmetry  axis  of  the  universe. 

As  the  polarimeter  sweeps  out  a  circle  of  constant  declination  night 
after  night,  we  can  add  the  output  records  synchronously  at  the  sidereal 
rate  and,  in  this  way,  perform  the  time  average  needed  to  overcome 
the  noise  temperature  of  the  receiver.   If  we  alternate  the  orientation 
of  the  radiometer  between  two  positions  U5  degrees  apart,  we  can  make 
the  two  orthogonal  measurements,  necessary  to  uniquely  specify  the 
polarization,  simultaneously. 

From  equation  2.l6  we  can  see  that  by  grouping  all  the  constants 

together,  the  result  can  be  expressed  more  simply  as 

T.  -  T0  =  (T  -  T  )     Ia  +   Bcos2(ar  -  t  +  6.,  )  +  C  cce  (a  -   t  +  6_)>   , 
1   2   ^  w   a'max  )         x  o      1J  v  o      2')      ' 

(2.17) 
with  relationships  between  the  coefficients  easily  derivable  from  the 
parent  equation.  This  form  displays  the  harmonic  dependence  upon  local 
sidereal  time  and  thus  the  dependence  upon  right  ascension.   In  general 
then,  each  measured  parameter  will  have  a  12  hour  and  a  2k   hour  harmonic 
component.   It  is  the  object  of  this  experiment  to  measure  the  polariza- 
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tion  parameters  as  well  as  possible  and  look  for  the  harmonic  dependence 
which  would  indicate  the  existence  of  a  large  scale  polarization  of 
cosmological  origin,  obeying  this  simple  model. 

2.3  THE  FARADAY  SWITCHED  POLAEIMETER 

The  device  used  to  perform  this  experiment  is  an  offspring  of  the 
Dicke  type  receiver  called  a  Faraday  switched,  polarimeter.   As  is 
shown  in  Figure  2.4,  instead  of  two  polarization  sensitive  antennas 
connected  to  the  receiver  through  a  single  pole  double-throw  switch, 
one  conical  type  polarization  insensitive  antenna  is  used.  This  horn 
feeds  into  an  axially  symmetric  chamber  containing  a  thin  ferrite  rod. 
Application  of  a  sinusoidal  magnetic  field  to  the  ferrite  causes  the 
incoming  radiation  to  be  rotated  back  and  forth  synchronously.   The 
chamber  is  followed  by  a  dual  mode  transducer  which  resolves  incoming 
radiation  into  two  orthogonal  components  and  thus  acts  as  an  analyzer. 
One  port  of  the  dual  mode  transducer  is  connected  to  the  receiver 
while  the  other  is  dumped  into  a  suitable  matched  load.  The  net 
effect  of  this  arrangement  is  the  same  as  was  obtained  using  t\>?o 
antennas,  only  now  some  of  the  matching  problems  have  been  eliminated. 
With  the  field  on  the  ferrite  in  one  direction,  radiation  from  one 
polarization  is  rotated  by  the  Faraday  effect  and  passes  through  the 
dual  mode  transducer  into  the  receiver.  If  the  field  is  then  reversed 
and  the  magnitude  has  been  adjusted  so  that  the  radiation  is  rotated 
90  degrees  in  going  from  one  state  to  the  other,  the  orthogonal 
polarization  will  now  be  presented  to  the  receiver.   In  this  way  a 
switched  signal,  proportional  to  the  temperature  difference  between 
the  tv/o  directions,  is  produced. 
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The  rest  of  this  chapter  will  be  devoted  to  this  device  and  its 
calibration.   Frequent  reference  will  be  made  to  Appendix  B  where 
the  detailed  specifications  and  measured  parameters  are  given. 

2.k     THE  ANTENNA 

A  conical  optimum  gain  horn  with  an  aperture  choke  for  suppression 
of  the  backlobes  was  used  in  this  experiment.*  Provision  was  made  for 
insertion  of  a  calibration  probe  in  the  side  of  the  horn  and  this 
feature  was  used  for  part  of  the  time.  Later  analysis  showed  that  the 
probe  increased  the  sensitivity  of  the  polarimeter  to  ground  radiation, 
an  effect  which  is  described  quantitatively  in  Appendix  Bl. 

The  full  width  of  the  main  lobe  between  half  power  points  was 
about  15  degrees,  giving  the  polarimeter  a  resolution  of  one  hour 
in  right  ascension. 

2.5  THE  SVHITCH 

The  heart  of  the  polarimeter  is  the  Faraday  switch  used  to  rotate 
between  orthogonal  polarizations.  Because  it  was  used  to  produce  the 
signal,  any  departure  from  ideal  operation  will  intimately  affect  the 
results  of  the  experiment.   In  this  section  I  will  present  those 
principles  of  operation  necessary  to  an  understanding  of  the  polarimeter, 
A  more  detailed  analysis  of  the  operating  principles  and  parameters  of 
the  switch  and  how  they  are  affected  by  outside  influences  is  given 
in  Appendix  B2. 

The  basic  principles  can  be  most  easily  understood  by  considering 
the  switch  as  an  ideal  Faraday  rotator.  As  shown  in  Figure  2.5,  the 
active  element  is  a  thin  cylindrical  piece  of  fcrrite  mounted  coaxially 
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in  a  section  of  circular  waveguide  around  which  is  wrapped  a  ceil  of 
wire.  Putting  a  current  through  the  coil  produces  a  longitudinal 
magnetic  field  in  the  ferrite  which  in  turn  produces  the  Faraday 
effect,  by  causing  the  magnetic  moments  of  the  ferromagnetic  ions  to 
align  with  the  field.  If  then,  a  circularly  polarized  wave  propagates 
along  the  axis,  its  velocity  is  either  decreased  or  increased  depending 
on  whether  its  helicity  has  the  same  or  opposite  sense,  relative  to 
the  precession  of  the  ferromagnetic  moments.  Since  a  plane  polarized 
wave  is  composed  of  equal  left  and  right  circularly  polarized  components, 
one  component  is  retarded  relative  to  the  other  and  the  plane  of 
polarization  is  rotated.  By  placing  an  oscillating  current  of 
sufficient  magnitude  on  the  driving  coil,  the  propagating  wave  can  be 
rotated  first  in  one  direction  by  k5   degrees  over  the  length  of  the 
ferrite,  and  when  the  current  is  reversed,  k^   degrees  in  the  opposite 
direction.  Now  if  the  rotator  is  followed  by  a  dual  mode  transducer 
which  acts  as  an  analyzer,  separating  the  radiation  into  two  orthogonal 
components,  one  obtains,  in  either  arm  of  the  transducer,  radiation 
which  is  modulated  at  the  switching  frequency.  The  amplitude  of  the 
modulated  signal  is  proportional  to  the  difference  between  two 
polarization  components  of  the  incoming  radiation,  referenced  to  a 
coordinate  system,  which  has  one  of  its  axes  bisecting  the  angle 
between  the  sensitive  directions  of  the  dual  mode  transducer.  The 
graphs  in  Figure  2.6  depict  the  result  for  incident  polarized  radiation 
in  two  situations.  In  one  case,  the  direction  of  polarization  is 
parallel  to  one  of  the  axes  of  the  reference  coordinate  system,  giving 
the  maximum  signal.   In  the  other  case,  the  incident  radiation  is  at 
an  angle  of  J+5  degrees  to  this  direction  or,  in  other  words,  parallel 
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to  the  direction  of  one  of  the  sensitive  directions  of  the  dual  mode 
transducer,  giving  zero  signal.  That  this  second  case  in  fact  gives 
zero  signal  comes  about  because  the  output  has  a  periodicity  at  twice 
the  switching  frequency.  Since  the  lock-in  extracts  only  the  first 
harmonic  of  the  switching  frequency,  this  contribution  is  zero. 

It  should  be  noted  that  I  am  looking  at  temperature,  a  quantity 
which  for  an  extended  blackbody  source  is  related  to  the  incident 
power  by 

P  =  kTAv,  (2.18) 

in  the  Rayleigh-Jeans  portion  of  the  spectrum.  Therefore,  I  am 
sensitive  to  a  quantity  which  is  proportional  to  the  square  of  the 
incident  field.   Because  of  the  lack  of  dependence  on  the  sign  of  the 
incident  field,  positions  of  the  polarimeter  which  differ  by  a 
rotation  of  180  degrees  about  the  ferrite  axis  are  equivalent,  while 
positions  90  degrees  apart  are  opposite  in  sign. 

In  practice  this  symmetry  is  not  perfect,  being  broken  by  an 
offset  due  to  absorption  and  emission  of  radiation  by  the  switch 
itself.  To  make  matters  worse, this  offset  is  a  funqtion  of  three 
major  perturbing  influences;  the  temperature  of  the  ferrite,  changes 
in  ferrite  magnitization  due  to  DC  fields  like  the  earth's  magnetic 
field,  and  mechanical  strains  transmitted  to  the  ferrite  through  its 
mountings.   In  order  to  reduce  the  effects  of  temperature,  the 
polarimeter  was  mounted  in  a  temperature  controlled  environment, 
where,  except  under  extreme  conditions  of  heating  or  cooling,  the 
temperature  was  maintained  to  about  one  degree  centigrade.  Effects 
from  magnetic  fields  were  removed  by  encasing  the  entire  microwave 
front  end  of  the  device  in  a  shield  constructed  of  extremely  high 
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permeability  material.   Finally,  mechanical  strains  were  kept  to  a 
minimum  by  mounting  the  axis  of  the  polarimeter  vertically,  and 
thus  parallel  to  the  gravitational  force. 

One  can  choose  to  take  the  signal  out  of  either  arm  of  the  dual 
mode  transducer,  since  these  outputs  give  the  same  signal,  differing 
only  by  180  degrees  in  phase  at  the  switching  frequency.   In  fact, 
it  would  be  possible  to  use  both  outputs  to  advantage  as  will  be 
discussed  later.  In  this  experiment  I  used  only  one  arm,  attaching 
a  small  standard  gain  horn,  aimed  at  the  sky,  to  the  other.   In  this 
way  the  unwanted  signal  was  discarded  and  the  unused  port  of  the 
switch  was  terminated  in  a  cold  load. 

Returning  to  the  subject  of  the  switch  offset,  though  I  have 
described  how  perturbing  influences  were  kept  from  altering  it,  I 
haven't  yet  explained  how  to  eliminate  it  as  an  instrumental  effect. 
The  easiest  way  to  accomplish  this  is  to  take  alternate  readings  of 
the  polarization,  separated  by  a  rotation  of  90  degrees  about  the 
ferrite  axis.  This  would  result  in  two  signals 


Sl  "  Tl  +  Toffset 


S2  "  "  Tl  +  Toffset' 


(2.19) 


and  by  subtracting  them  we  would  obtain 

S1-  S2  =  2Tr  (2.20) 

a  quantity  in  which  the  offset  has  been  subtracted  out.   Of  course, 
this  procedure  depends  on  being  able  to  shield  out  the  effects  of  the 
earth's  magnetic  field  and  the  perturbiag  influence  of  gravitational 
strains.  Having  the  polarimeter  axis  vertical  keeps  the  gravitational 
forces  constant  with  rotation. 
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2.6  THE  CIRCULATOR  AND  ISOLATOR 

These  components  are  designed  to  provide  isolation  between  the 
microwave  mixer  and  the  back  of  the  switch.  Because  the  local 
oscillator  is  not  perfectly  monochromatic  and  still  has  a  power 
contribution  at  the  sideband  frequencies,  isolation  is  made  necessary 
to  prevent  power  leakage  back  up  the  waveguide.   For  the  Gunn 
oscillator  used,  the  power  in  the  sidebands  was  on  the  order  of 
-1^0  db  below  the  center  frequency  power.   As  this  is  a  factory  claim 
and  wasn't  directly  measured,  adding  a  factor  of  ten  might  be 
appropriate;  so  -130  db  might  be  more  comfortable.   For  input  powers 

of  between  one  and  two  millawatts  this  means  that  there  may  be  as 

-1 6 

much  as  2  x  10    watts  left  in  the  sidebands .   It  might  not  seem  like 

much  but  being  in  the  Rayleigh-Jeans  portion  of  the  spectrum  and  with 
a  bandwidth  of  9&  MHZ,  we  find  that  a  millidegree  corresponds  to 

-1  Pi 

l.k   x  10    watts.   In  short,  we  may  find  as  much  as  100  rnillidegrees 
worth  of  power  leaking  out  of  the  mixer  in  the  sidebands  we  are 
sensitive  to. 

Ordinarily  this  would  probably  be  no  problem,  but  the  reflection 
coefficient  of  the  ba  ck  of  the  switch  is  a  strong  function  of  magneti- 
zation.  For  that  reason  one  doesn't  want  to  provide  it  with  a  nice 
source  of  power  at  the  sideband  frequencies,  which  it  could  modulate 
through  reflection  and  send  back,  masquerading  as  real  signal. 

To  provide  isolation,  a  circulator  and  isolator  were  inserted 
between  the  switch  and  mixer.  Figure  2.7  diagrams  the  operating 
characteristics  of  the  circulator.  Power  leaving  the  switch  is  passed 
to  the  mixer.  Power  from  the  mixer  is  passed  to  the  third  port  which 
is  terminated  in  a  standard  gain  horn  aimed  at  the  sky,  and  a  cold 
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5  K  from  the  sky  illuminates  the  back  of  the  switch,  giving  it  little 
signal  to  modulate  by  reflection.   About  30  db  worth  of  isolation  is 
provided  in  this  way. 

Another  kO   db  is  provided  by  the  isolator,  which  is  inserted 
between  the  circulator  and  the  mixer.  It,  in  fact,  is  also  a  circulator 
but  with  its  third  port  terminated  internally  in  a  matched  load.  The 

addition  of  these  two  devices  reduces  the  residual  local  oscillator 

-23  -5 

power  reaching  the  switch  to  about  10    watts  or  10   mdeg.  K, 

a  much  more  acceptable  level. 

Because  these  two  devices  depend  on  ferrites  for  their  operation, 

it  was  thought  that  perhaps,  like  the  switch,  they  would  be  subject 

to  variations  in  their  operation  due  to  stray  magnetic  fields.   It 

might  not  have  been  as  severe  in  this  case,  as  neither  of  these 

devices  uses  an  oscillating  field  to  produce  a  signal.  None  the  less, 

they  pass  the  offset  produced  by  the  switch  and  if  their  absorption  is 

changed  by  a  large  enough  amount,  they  could  produce  an  orientation 

dependent  signal.   For  this  reason  they  were  tested  by  an  application 

of  field  along  the  axes  of  their  ferrite  devices.  This  should  be 

their  most  sensitive  directions.  No  measurable  effect  was  observed 

for  fields  up  to  11  Gauss. 
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2.7  THE  RECEIVER 

A  standard  off  the  shelf  model  heterodyne  receiver,  having  a 
separate  double  balance  mixer  and  IF  amplifier  was  used  at  9»37  GHz. 
The  width  of  each  sideband  was  U8  MHz  ,  giving  a  total  receiver  band- 
width of  96  MHz.   Local  oscillator  power  was  provided  by  a  Gunn 
oscillator  coupled  through  a  frequency  meter  and  attenuation  pad. 

The  only  unusual  operating  characteristic  of  this  receiver  was 
a  bad  impedance  match  between  the  mixer  preamplifier  and  IF  amplifier. 
If  a  resonant  length  of  cable  were  used  and  the  system  perturbed  by 
some  electrical  interference,  oscillatory  breakdown  could  occur. 
By  judicious  choice  of  cable  lengths  and  careful  shielding  against 
interference,  the  problem  was  eliminated.   Because  a  chart  record  of 
the  lock-in  output  was  kept  and  because  this  problem  gave  a  striking 
signature,  those  times  that  it  did  occur  (as  the  result  of  extreme 
electrical  interference  in  thunderstorms,  for  example)  were  easily 
seen  and  the  affected  data  eliminated. 

One  modification  was  made  to  the  standard  IF  amplifier  and  that 
was  the  addition  of  a  ENC  terminal  to  allow  measuring  the  DC  voltage 
at  the  second  detector  before  it  was  capacitively  coupled  to  the 
video  stage.  As  will  be  explained  in  the  next  section,  this  was 
needed  for  the  calibration  procedure  used. 
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2.6  CALIBRATION 

Because  there  were  no  standard  polarized  blackbody  sources 
available  for  calibration  at  3.2  cm,  another  method  had  to  be  employed. 
The  procedure  used  involves  two  steps.   First,  the  change  in  voltage 
at  the  second  detector  as  a  function  of  antenna  temperature  is  deter- 
mined, followed  by  a  measurement  of  the  gain  through  the  rest  of 
the  system.  ' 

The  calibration  was  begun  by  measuring  the  DC  voltage  at  the 

second  detector  with  the  polarimeter  looking  at  the  sky.  The  size  of 

this  voltage  is  proportional  to  the  power  entering  the  horn  from  the 

sky  plus  some  offset  due  to  receiver  noise  and  biasing  of  the 

electronics.  Next,  a  piece  of  eccosorb  microwave  absorber,  a  good 

unpolarized  blackbody  source  at  ambient  temperatures,  was  placed  over 

the  antenna  and  the  voltage  again  read.  This  time  the  voltage  is 

proportional  to  the  power  from  the  eccosorb  plus  the  same  offsets 

that  were  present  before.  In  this  way  we  have  measured  the  following 

two  quantities: 

v  i   =  C,  (T  .   +  T  ._   ,  ) 
sky    1  v  sky    offset7 

V      v  -  C,  (T      .  +  T  _-  .),  (2.21) 

eccosorb    1  x  eccosorb    offset''  v    ' 

where  C  is  a  constant  dependent  on  the  gain  of  the  receiver.  Combina- 
tion of  these  two  allows  one  to  determine  C  and  thus  the  response 
in  Volts/  K  of  the  receiver  to  this  point.  Now  imagine  that  the 
polarimeter  is  looking  at  an  extended  blackbody  source  which  is  100$ 
polarized,  such  that  the  temperature  in  the  direction  of  polarization 
is  equal  to  the  temperature  of  the  eccosorb  and  such  that  the 
temperature  in  the  orthogonal  direction  is  equal  to  the  temperature 


35 


of  the  sky.   If  this  source  were  aligned  with  the  symmetry  axis  of  the 

switch,  it  would  generate  a  signal  at  the  second  detector  at  the 

switching  frequency,  like  that  in  Figure  2.8.   In  this  way  one  can  say 

that  a  polarized  source  of  magnitude  T       ,  -T  ,   will  produce  a 
*  eccosorb  sky     L 

signal  at  the  second  detector  whose  effective  first  harmonic  peak-to- 
peak  value  is  l.lGk   (V      .-V  ,  ),  where  the  calibration  value 

eccosorb  sky 

derived  in  Appendix  B2  has  been  used.   Redefining  C  to  accommodate 

the  calibration  factor,  it  can  be  written 

1.164  (V  ,  -V      .  ) 
n sky  eccosorb7  ,    s 

sky  eccosorb 

where  a  typical  value  for  V  „  -V      ,  was  7.0  mv  measured  to  an 

sky  eccosorb 

accuracy  of  .1  mv  and  for  T  ,  -T      ,  was  -2Q0.0  K  measured  to 

sky  eccosoro 

an  accuracy  of  .25  K. 

Next,  I  insert  a  calibration  signal  at  the  second  detector.  This 
is  a  sinewave  of  known  peak-to-peak  amplitude,  which  is  derived  from 
the  audio  oscillator  through  a  special  circuit.  With  the  lock-in 
phase  adjusted  for  maximum  signal,  the  output  voltage  is  read  and  a 


second  parameter  C„,  which  gives  the  gain  of  the  rest  of  the  receiver 


is  computed: 

C2  =  ^liR  .  (2.23) 

calibration  signal 

The  product  of  C?  and  C.  then  gives  the  gain  of  the  entire  polarimeter. 

In  practice, the  measurement  was  done  utilizing  the  output  of  the 

entire  data  collection  system  rather  than  just  that  of  the  lock-in. 

In  this  way  the  overall  gain  of  the  entire  system  could  be  measured 

and  individual  intermediate  steps,  where  human  reading  errors  could 

multiply,  were  minimized. 
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Although  the  temperature  of  the  eccosorb  was  always  known,  that  of 
the  sky  couldn't  be  with  great  certainty.   For  this  reason,  I  picked  a 
standard  temperature  of  5  K  for  the  sky  for  purposes  of  calibration, 
even  though,  typically,  the  actual  temperature  will  fluctuate  from 
5-10  K.   It  was  felt  that  the  error  introduced  by  this  assumption 
was  only  on  the  order  of  5/300  or  about  2cj0,   which,  compared  to  the 
uncertainty  due  to  receiver  noise  alone,  is  negligible. 

After  a  few  calibrations  it  became  clear  that  the  overall  gain 
stability  of  the  system  was  good  enough  thac  a  calibration  on  the 
order  of  once  a  week  was  enough  to  ensure  good  accuracy.  Table  2.1 
is  a  list  of  calibrations  taken  during  the  data  collection  period. 
The  units  in  which  C  C?  is  expressed  are  counts  per  second  at  the 
output  of  the  voltage-to-frequency  converter  per  millidegree  Kelvin 
at  the  polarimeter  input.   The  output  of  this  device  is  directly 
proportional  to  the  lock-in  output  and  was  used  to  make  integration 
of  the  signal  and  digitization  of  the  result  much  easier.  The  entries 
are  segregated  into  two  groups  because  work  on  the  receiver  changed 
the  overall  gain  midway  through  the  data  taking  period.   Figure  2.9 
plots  the  distributions  of  C.  C?  for  the  two  groups,  which  appear 
approximately  normally  distributed.   It  is  interesting  to  note  that 
if  the  data  had  been  interpreted  using  the  mean  of  each  group  of 
calibrations,  the  expected  error  in  the  calibration,  equal  to  the 
standard  deviation  of  the  measurements  divided  by  the  mean,  would  be 
about  2.8$. 

As  a  final  comment,  I  would  like  to  discuss  the  overall  effective- 
ness of  this  method  considering  that  there  was  no  true  polarized 
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TABLE  2.1 
POLARIMETER  CALIBRATIONS 


First  Group 

Date 

C  C2    (CNTS/mdeg  K) 

Date     . 

C-Cg   (CNTS/mdeg  K) 

8/11/72 

1+0.50 

11/20/72 

40.31 

8/2U/72 

39.91 

II/29/72 

40.91 

9/7/72 

42.28 

12/7/72 

39.81 

9/18/72 

40.75 

12/14/72 

39.84 

9/26/72 

40.57 

12/29/72 

39-08 

10/2/72 

4l.79 

1/V73 

40.43 

10/10/72 

42.85 

1/7/73 

39.11 

10/23/72 

4l.  96 

l/n/73 

41.64 

10/30/72 

39-82 

1/18/73 

42.61 

11/9/72 

41.55 

1/26/73 

41.75 

OjCg  =  40.83 

orc1c2  =  1.14 

a          1  — 

2  =    .028 

Second  Group 


Date 

C   Cg    (CNTS/mdeg  K) 

Date 

C-Cg   (CM'S/mdeg  K) 

2/4/73 

43.93 

4/24/73 

45.00 

2/13/73 

44.42 

5/6/73 

44.49 

2/20/73 

42.40 

5/29/73 

43.32 

2/28/73 

42.72 

6/11/73 

44.81 

3/9/73 

44.12 

7/4/73 

45.63 

3/20/73 

44.21 

7/20/73 

44.07 

3/28/73 

44.16 

8/1/73 

45.24 

4/6/73 

44.48 

8/10/73 

47.47 

4/16/73 

43.00 

C1C2  =  UJ4*32 

CTc1c2  =1.19 

V/c. 

c„  =  0.27 

1    <= 
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standard  against  which  to  calibrate.  As  a  means  to  justify  the  method, 
I  made  a  comparison  between  the  results  given  by  my  polarimeter  and  the 
results  given  by  a  polarization  sensitive  power  radiometer,  both 
measuring  the  same  polarized  source.  Used  as  a  source  was  a  klystron 
radiating  through  a  rectangular  standard  gain  horn  suspended  four 
stories  above  the  measuring  devices  and  giving  an  almost  fully 
polarized  signal.   A  reading  with  the  polarimeter  was  sandwiched  in 
between  two  readings  taken  with  the  power  radiometer.  Agreement 
between  the  two  radiometer  readings  was  2%,   while  agreement  between 
their  mean  and  the  polarimeter  reading  was  U-%.     Because  the  power 
radiometer  and  the  polarimeter  used  most  of  the  same  components  (the 
only  difference  being  the  addition  of  a  second  dual  mode  transducer 
and  comparison  load  between  the  antenna  and  the  mouth  of  the  switch), 
changing  configuration  meant  partial  dismantling  of  the  apparatus. 
During  this  time  great  care  was  taken  to  maintain  alignment  with  the 
source,  but  not  all  the  microwave  plumbing  was  uniform.   For  this 
reason  it  was  felt  that  the  errors  involved  were  ones  of  alignment 
and  not  ones  of  calibration.   The  error  in  each  measurement  due  to 
receiver  noise  was  only  about  15  mdeg.K,  while  the  signal  was  about 
280  K.   Given  this,  the  discrepancy  between  the  two  power  radiometer 
measurements  must  also  have  been  due  to  alignment,  as  the  calibration 
constant  only  varied  about  .1%   between  the  two  runs. 
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III .   PROCEDURE 

3.1  EXPERIMENTAL  PROCEDURE 

For  the  reasons  explained  in  previous  sections,  data  were  taken 
by  the  polarimeter  with  its  axis  mounted  vertically  and  observing  the 
zenith.  Rotation  of  the  earth  sweeps  out  a  circle  of  constant 
declination  in  the  sky.   By  rotating  the  polarimeter  by  k-5   degrees 
between  observations,  measurement  of  the  Q  and  U  Stokes  parameters 
can  be  made. 

'There  is  one  additional  source  of  systematic  error  inherent  in 
this  method.  This  source  can  be  eliminated  by  data  taking  procedures. 
In  Chapter  2  and  Appendix  B2,  I  discuss   at  length  the  problems  of 
external  parameters  influencing  the  offset  of  the  Faraday  rotation 
switch.  Aside  from  these  problems,  it  is  easy  to  show  that  any 
offset  at  all  is  harmful.   For  example,  for  a  switch  offset  of  50 
millidegrees  entering  the  receiver,  any  gain  fluctuations  on  the 
order  of  2%  will  produce  systematic  effects  of  around  a  millidegree 
in  the  data,  an  unacceptable  level.  To  remove  this  we  must  find  some 
way  of  eliminating  the  offset  from  the  data  in  such  a  way  that  the 
effects  of  nonrandom  gain  fluctuations  will  be  eliminated. 

A  possible  solution  utilizes  the  symmetry  properties  of  the 

switch  itself.  With  the  polarimeter  aligned  with  one  polarization 

axis  north-south,  it  produces  a  signal  proportional  to 

Tmo-T  +T   „  ,  .  Rotating  the  axis  of  the  polarimeter  by  ninety  degrees 
No     ejY!     oil  set 

produces  a  signal  T„TT-  !„_+  T  „_   ,  where  the  contribution  due  to 
*  EW   NS   offset 

T  -.„  ,  doesn't  change.   It  is  a  systematic  effect  in  the  -oolarimeter 
offset  °  * 

itself,  and  is  independent  of  external  orientation  if  the  switch  is 
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properly  shielded.  The  readings  taken  in  the  two  positions  are 
subtracted  to  give 

S  -S 

— —  =  T   -  T 

2        NS   EW  ' 

from  which  the  offset  has  been  eliminated.   Of 'course,  how  well  this 

subtraction  works  depends  on  how  constant  the  offset  is  over  the 

rotation  timescale.  Keeping  the  integration ,times  short  reduces 

sensitivity  to  longer  term  drifts  in  the  offset. 

Combining  the  idea  of  rotating  the  polarimeter  by  ninety  degrees 
between  equivalent  readings  with  the  need  to  take  readings  separated 
by  forty-five  degrees,  one  arrives  at  a  natural  solution  to  both 
problems;  that  is,  by  rotating  by  forty-five  degree  steps  around  the 
compass.   In  this  way  one  measures  both  Stokes  parameters  and  for 
each  measurement  there  are  two  others  separated  from  it  by  ninety  degrees, 
to  subtract  out  switch  offset. 

The  polarimeter  was  mounted  on  a  rotating  platform  with  magnetically 
actuated  reed  switches  located  at  forty-five  degree  intervals  around 
the  periphery,  and  labeled  one  through  eight.  Data  taking  was  con- 
ducted by  rotating  the  table  in  a  clockwise  direction  until  it  actuated 
a  switch,  stopping  and  integrating  for  five  minutes  at  that  position, 
and  then  recording  the  result  and  moving  on  to  the  next  position. 
When  the  recording  at  position  eight  was  completed,  the  table  would 
reverse  direction  and  return  to  position  one. 

To  obtain  the  most  useful  results,  it  was  felt  that  the  following 
experimental  quantities  had  to  be  measured: 
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1.  The  output  of  the  integrator,  which  sums  the  polarimeter 
output  for  five  minutes „ 

2.  The  time  of  the  observation.   This  quantity  can  be  easily 
converted  to  local  mean  sidereal  time  which  is 
identical  to  the  right  ascension  of  the  observation. 

3.  The  angular  position  of  the  polarimeter,  which  determines 
which  Stokes  parameter  one  is  measuring  and  with  which 
sign. 

k.      The  temperatures  of  the  switch,  the  inside  of  the  observa- 
tion shed,  and  the  outside  environment.  The  first 
of  these  was  needed  to  correlate  the  switch  offset 
with  temperature.   I  also  kept  track  of  the  inside 
shed  temperature  to  insure  that  my  temperature 
controller  was  working  properly.  In  the  event  that- 
ground  radiation  became  a  large  problem,  measuring 
the  temperature  of  the  surrounding  buildings  would 
help  in  removing  the  effect. 
The  apparatus  involved  in  taking  and  recording  these  observations 
is  the  subject  of  the  next  few  sections.   Figure  3.1  is  a  diagram  of 
the  entire  system.   It  was  felt  that  by  automating  as  much  as 
possible,  the  chances  for  human  error  could  be  reduced,  especially 
in  an  undertaking  of  this  type  where  the  observations  extend  over 
many  months. 
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3.2  PREPARATION  OF  SIGNAL 

The  output  of  the  lock-in  amplifier  was  recorded  in  two  ways. 
First  the  output  itself  was  recorded  directly  so  that  any  abnormal- 
ities in  function  or  signal  received  could  be  detected  and  evaluated. 
Two  examples  which  were  easiest  to  detect  in  this  way  were  rain  and 
local  thunderstorm  activity.  The  most  dependable  way  found  for 
recording  was  a  hand  wound, galvanometer  movement,  Esterline  Angus 
chart  recorder  which  could  be  left  unattended  for  two  to  three  days 
at  a  time.  Combined  with  a  marking  device  which  delineated  the 
integration  periods  in  each  position  and  the  hours  of  the  day,  I  was 
able  to  relate  a  disturbance  in  the  apparatus  to  a  recorded  data 
point . 

To  integrate  the  output  of  the  lock-in,  I  used  a  prescaler  and 
voltage -to -frequency  converter  which  changed  a  -10  to  +10  volt 
output  to  a  0  to  100,000  Hz  signal.  Integration  was  performed  by 
a  remotely  controlled  scaler,  which  was  capable  of  integrating  the 
signal  for  periods  from  zero  to  9999  seconds,  could  divide  the  input 
by  powers  of  ten,  and  could  display  six  digits  of  output  information 
in  binary  coded  decimal  form.   This  result  was  presented  in  parallel 
form  to  a  parallel  to  serial  converter  which  prepared  it  for 
recording. 
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3.3  TIMING 

Time  information  was  provided  by  a  clock  built  by  E.  Groth  and 
used  primarily  for  pulsar  timing.  Needless  to  say  its  micro-second 
accuracy  was  not  needed,  but  it  did  provide  time  in  days,  hours, 
minutes  and  seconds  of  Eastern  Standard  Time  relative  to  epoch 
1970.0.  Upon  completion  of  each  integration  period,  the  output 
of  the  clock  was  stored  and  the  result  presented  to  the  parallel- 
to-serial  converter  for  recording.  This  provided  a  convenient 
and  reproducable  point  from  which  to  calculate  the  time  of  the 
middle  of  each  integration  run. 

3.k     TEMPERA r ORE  MEASUREMENT 

The  temperature  at  key  locations  was  measured  using  four 
linearized  thermistor  composites  connected  to  a  measuring  circuit 
and  digital  voltmeter  through  a  multiplexing  scanner.  During  each 
run  one  of  the  thermistor  probes  was  measured  and  recorded.  The 
output  of  the  digital  voltmeter  and  the  position  of  the  scanner  were 
presented  to  the  parallel -to-serial  converter  for  recording.   It  was 
thought  that  this  information  might  be  used  to  correlate  out  the 
temperature  dependent  effects . 
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3.5  RECORDING  OF  DATA 
The  data  outputs  are  ultimately  recorded  in  IBM  compatible  form 

on  digital  tape.   The  intermediary  in  the  process  is  the  parallel-to- 
serial  converter.  When  strobed,  it  takes  all  the  data  presented  to  it 
in  parallel  form  and  feeds  it  serially  one  byte  at  a  time  to  an 
incremental  nine  track  800  byte  per  Inch  tape  unit.   In  addition,  it 
controls  the  blocking  of  the  data  into  files  'and  records  to  facilitate 
analysis.  This  device  plus  the  scaler  and  dvm  scanner  were  constructed 
by  Karl  C.  Davis. 

3.6  EXPERIMENT  CONTROL 
Overall  control  and  sequencing  of  the  experiment  was  performed  by 

an  experiment  control  box,  the  function  of  which  was  to  command  the 
polarimeter,  the  integrator  and  parallel-to-serial  converter  and 
sequence  their  functions  properly.   The  order  of  operations  while  the 
experiment  was  being  conducted  was : 

1.  Rotate  to  a  new  position.   If  at  position  eight,  rewind 
until  reaching  position  eight  again  and  then  step 
forward  to  position  one . 

2.  Start  integration. 

3.  Stop  integration  and  record  data. 
k.     Recording  finished,  rotate  to  a  new  position. 

Cycle  time  for  the  apparatus  with  a  three  hundred  second  integration 
time  was  three  hundred  and  three  seconds  except  when  rewinding,  in  which 
case  it  was  three  hundred  and  twenty  seconds.   The  device  would  maintain 
this  timing  without  losing  a  second  for  a  week  or  more,  making  it  very 
easy  to  spot  malfunctions  and  timing  errors. 
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All  switching  in  the  system  is  either  low  voltage  or  solid  state 
to  avoid  radio  frequency  interference  with  the  sensitive  receiver. 
For  this  reason.  AC  switching  for  the  polarimeter  rotation  motors  is 
tied  to  the  zero  crossing  of  the  line  voltage,  using  a  semi- 
conductor triac. 
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IV.   ANALYSIS  AND  INTERPRETATION 

k.l     DATA  EDITING 

Data  was  taken  for  this  experiment  from  1  February  1972  until 
10  August  1973.  The  data  before  9  August  1972  was  mainly  for  test  and 
evaluation  of  the  device  and  has  been  discarded  out  of  hand.  The  data 
since  that  time  is  considered  suitable  for  analysis  and  has  been 
searched  for  the  effects  of  cosmological  polarization.  The  twelve 
month  data  set  contains  approximately  9^j000  individual  polarization 
measurements,  of  which  about  half  were  devoted  to  each  Stokes 
parameter.  As  in  any  experiment  of  this  length,  bad  weather,  equip- 
ment malfunctions,  and  human  errors  took  their  toll,  resulting  in 
discarded  readings. 

My  guiding  principle  in  throwing  away  data  was  not  to,  if  at  all 
possible.   When  necessary  it  would  have  been  best  to  completely  auto- 
mate the  procedure  in  order  to  avoid  human  intervention  with  its 
attendant  bias.  But  at  certain  points,  in  particular  the  screening 
of  chart  recordings  of  the  lock-in  output,  there  was.  no  other  way. 
In  general  there  were  four  sources  of  abnormality  in  the  lock-in 
output : 

1.   The  most  common  and  troublesome  was  rain.  The  polarimeter 

was  insensitive  to  clouds  and  other  forms  of  bad  weather, 
but  precipitation  would  collect  on  the  mylar  window  covering 
the  apparatus  and  appear  as  a  strongly  polarized  signal. 
This  would  many  times  drive  the  lock-in  off  scale. 
Data  showing  this  type  of  behavior  was  rejected  along 
with  enough  data  on  either  side  to  insure  that  all 
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contamination  was  removed. 

2.  The  next  most  frequent  source  of  interference  was  thunder- 

storm activity,  especially  in  the  summer  months.  This 
was  characterized  by  strong  and  repeated  noise  spikes 
in  the  output  signal,  which  caused'  the  data  to  be 
discarded.   Isolated  noise  spikes  were  not  thrown 
away  by  manual  means.   Later  stages  of  automated 
editing  evaluated  their  effect  and  provided  more 
impartial  criteria  for  their  elimination. 

3.  During  the  summer  months  when  the  sun  reaches  altitudes 

close  to  the  zenith,  its  presence  as  a  hot  source 
produces  a  lock-in  signal  when  acting  on  the  non- 
symmetric  antenna  gain  pattern.  Again  this  data 
was  not  deleted  manually  but  was  handled  later  when 
all  data  with  the  sun  above  the  ground  shield  and 
geometrically  illuminating  the  antenna  were  discarded 
k.      Finally,  there  would  be  an   occasional  malfunction  in  the 
receiver  and  its  associated  equipment,  which  would 
usually  be  indicated  by  a  lack  of  receiver  noise  in 
the  output  trace.  This  was  easily  noticed,  the 
receiver  repaired  and  the  data  removed. 
A  miscellaneous  problem  that  does  not  fit  into  any  of  the  above 
categories,  but  which  at  one  point  had  me  convinced  that  the  apparatus 
was  broken,  was  that  of  a  small  very  inconspicuous  spider.   By  slowly 
crawling  around  in  the  antenna  throat  he  did  a  great  imitation  of 
serious  receiver  gain  drifts.  Naturally  his  virtuoso  performance 
was  eliminated. 
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Removal  of  bad  data  constituted  the  first  stage  of  editing.  Next 
I  searched  for  timing  errors  to  make  sure  that  the  automated  data 
collection  apparatus  was  functioning  properly  and  that  the  integration 
periods  were  uniform.  These  events  fell  into  two  categories,  clock 
errors  and  actual  erratic  operation.   In  the  former  case  the  equipment 
would  be  operating  normally,  but  the  clock  would  indicate  the  wrong 
time.   Resetting  of  the  clock  for  pulsar  timing  was  the  most  common 
cause  of  this  and  would  appear  as  a  step  function  in  the  cycle, 
which  would  be  followed  by  a  step  with  the  opposite  sign  when  the 
correction  was  completed.  The  cycle  time  of  the  apparatus  was 
consistent  enough  to  enable  one  to  reconstruct  the  correct  time 
even  when  the  step  remained  for  several  hours. 

Harder  to  correct  was  difficulty  due  to  erratic  behavior  of  the 
apparatus,  which  was  usually  associated  with  thunderstorm  activity 
or  other  forms  of  electrical  interference.  In  most  of  these  cases, 
where  the  orderly  sequence  of  readings  could  not  be  restored,  the 
errors  were  left  for  a  later  stage  where  they  would  automatically 
be  eliminated. 

On  the  average  of  every  one  or  two  weeks, the  cycle  would  lose  a 
second  or  two  in  phase  for  no  apparent  reason.  No  attempt  was  made 
to  correct  for  this.   In  later  processing,  isolated  errors  in  timing 
of  up  to  five  seconds  were  allowed  to  pass.  A  change  in  the 
integration  period  of  up  to  five  seconds  means  an  error  in  the 
measurement  of  under  two  percent  or  far  less  than  the  uncertainty 
due  to  receiver  noise  alone. 
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As  a  final  editing  procedure,  the  problem  of  isolated  interference 
spikes  in  the  data  was  considered.   Preliminary  examination  showed  that 
after  removal  of  switch  offset,  the  averaged  signal  would  be  within 
10  millidegrees  of  zero.  The  minimum  detectable  signal  or  standard 
deviation  of  the  noise  for  a  single  processed  reading  was  about  15 
millidegrees  based  on  measurements  of  scatter  in  the  data.   The 
clipping  level  for  rejection  of  spikes  was  set  three  full  standard 
deviations  above  the  10  millidegree  level,  at  56  millidegrees  Kelvin. 
Only  387  points  had  to  be  eliminated  in  this  way. 

In  all,  out  of  about  9^4,000  points  collected,  only  38,029  points 
or  about  kO.k'fo   of  the  data  were  actually  used.  Unfortunately  they 
were  not  uniformly  distributed  throughout  the  year  for  two  reasons. 
First,  more  data  get  thrown  out  in  summer  when  the  sun  is  higher  in 
the  sky.  Worsening  the  effect  was  a  very  wet  spring  and  early  summer 
in  the  Princeton  area,  combined  with  unstable  air  and  much  electrical 
storm  activity. 
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k.2     PRELIMINARY  PROCESSING 

The  final  stage  of  data  preparation  is  carrying  out  the  subtraction 
(described  under  experimental  procedure)  for  removal  of  the  switch 
offset.  At  this  point  we  have  two  roughly  equal  sets  of  measurements 
of  the  two  Stokes  parameters  describing  the  polarization.  Each  of 
these  sets  contains  measurements  at  four  of  the  polarimeter  positions. 
One  set  comprises  all  the  odd  numbered  switch  positions  and  the  other 
has  the  even,  as  in  Figure  k.l.     For  the  Q  parameters,  that  position 
closest  to  North  or  position  two,  was  called  positive  in  accordance 
with  normal  convention.   It  was  found  by  checking  with  a  polarized 
source  placed  due  North  of  the  apparatus  that  this  position  actually 
pointed  11.75  degrees  to  the  Northeast  of  North.  I  will  still  call 
this  the  Q'  parameter  (primed  to  distinguish  it)  and  correct  the 
discrepancy  at  a  later  stage.   For  the  U'  parameters,  the  positive 
direction  is  that  closest  to  the  Northeast  or  position  three. 
Positions  180  degrees  relative  to  two  and  three  have  the  same  sign 
by  symmetry,  while  those  90  degrees  away  have  the  opposite  sign. 

To  perform  the  subtraction,  each  reading  was  assigned  its 
characteristic  sign  and  then  adjacent  pairs  of  measurements  within 
each  set  were  averaged  as  in  Figure  k.2.      The  time  of  observation 
assigned  to  each  averaged  pair  is  the  average  of  the  midpoints  of 
each  integration  period. 

If  we  look  at  this  operation  from  another  point  of  view,  we  see 
that  what  we  are  doing  is  multiplying  the  data  by  a  periodic  signal 
of  frequency  twice  that  of  the  polarimeter  rotation  frequency.  This 
is  exactly  what  is  done  in  the  mixer  stage  of  a  lock-in  amplifier. 
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POLARIMETER  INTEGRATOR 
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Because  we  do  the  Q'  and  U'  parameters  separately,  we  extract  both 
phases  of  the  signal  simultaneously  just  as  in  a  two  phase  or  vector 
lock-in.  To  visualize  things  a  little  more  clearly,  consider  the 
effect  on  incident  DC  polarized  radiation  as  the  polarimeter  is 
rotated.  The  rotation  causes  a  sinusoidal  variation  in  intensity. 
Because  of  the  180  degree  symmetry  of  the  switch,  this  signal  has  twice 
the  rotation  frequency  of  the  polarimeter.  Now  by  multiplying  by  a 
periodic  signal  the  way  we  have  done,  we  demodulate  the  signal  and 
extract  the  fundamental  Fourier  component  for  each  phase,  exactly  as 
is  done  in  a  lock -in.  A  double  lock-in  technique  of  this  type  gives 
good  immunity  against  effects  from  gain  fluctuations  and  drifts  in 
the  switch  offset  which  are  on  a  time  scale  longer  than  half  the 
polarimeter  rotation  period,  or  twenty  minutes. 

1+.3  ANALYSIS 

I  now  have  two  sets  of  independent  measurements,  the  Q,'  and  U' 
parameters,  each  with  a  time  resolution  of  about  ten  minutes  in 
sidereal  time.   One  might  think  that  this  could  be  interpreted  as  a 
spatial  resolution  in  terms  of  right  ascension.  This  is  not  true 
because  of  the  low  resolution  of  the  antenna,  which  is  only  about  one 
hour  in  right  ascension.   To  bin  the  data  on  any  sidereal  time  scale 
shorter  than  this  is  to  attempt  to  display  resolution  which  does  not 
exist.  The  next  logical  step  is  to  bin  the  data  at  its  true  resolution 
limit,  one  hour  in  right  ascension  (local  sidereal  time). 

In  order  to  check  for  serious  solar  effects,  the  binning  was  also 
carried  out  in  solar  time.  After  binning,  Q'  and  U'  were  transformed 
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to  the  true  Q,  and  U  parameters  before  further  analysis.  The  results 
of  the  binning  are  shown  in  Figure  k.3   for  solar  time  and  Figure  h.k 
for  sidereal  time.  The  temporal  or  spatial  center  of  each  bin  may  be 
obtained  by  subtracting  thirty  minutes  from  the  bin  number.  In  other 
words,  bin  four  contains  all  the  data  between  three  and  four  hours. 

In  each  set  of  data  there  is  a  DC  level,  which  when  combined  from 
both  sets  gives  a  resultant  DC  polarization  yector  with  magnitude  of 
1.10  ±  .lk   millidegrees  Kelvin  and  direction  116.3*+  ±  3.6U  degrees  from 
North.   On  the  reciprocal  of  this  bearing  is  the  most  prominent  feature 
on  the  horizon,  a  tower  containing  elevator  motors  on  the  roof  of  Jadwin 
Physics  Laboratory.  When  reviewing  the  data  in  smaller  segments,  it  was 
found  that  this  level  wasn't  constant  throughout  the  year  of  observation. 
During  the  first  two  months,  the  DC  vector  was  6.6  ±  .28  millidegrees  with 
a  direction  of  150.U6  ±1.21  degrees.   Other  than  this,  no  other  abnormal- 
ities were  found.  Closer  examination  showed  that  the  change  took  place 
sometime  close  to  the  tenth  of  October.  Checking  my  notebook,  I  dis- 
covered that  on  the  fifteenth  of  October  I  had  decided  that  the  calibra- 
tion probe  installed  in  the  antenna  wasn't  being  used  and  should  be  removed. 
It  appeared  that  the  presence  of  a  stub  in  the  horn  had  perturbed  the 
gain  pattern  enough  to  significantly  alter  the  response  to  ground 
radiation,  in  effect,  making  the  polarimeter  six  times  more  sensitive 
to  the  elevator  tower  radiation.   The  discrepancy  in  angle  of  about 
+36  degrees  can  be  traced  to  the  fact  that  the  calibration  probe  was 
located  approximately  -kO   degrees  from  the  positive  symmetry  plane  of 
the  switch.  When  the  polarimeter  was  aligned  at  150  degrees,  the 
probe  was  aimed  at  the  tower.  A  check  was  performed  with  a  slotted 
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3.2   POLARIZATION  DATA 

FOLDED  INTO    SOLAR  BINS 

BIN  WIDTH  =  1hr 
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3.2    POLARIZATION  DATA 
FOLDED    INTO   SIDEREAL  BINS 
BIN   WIDTH  =  1hr 
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plate  providing  a  polarized,  source  to  determine  whether  or  not  the  probe 
changed  the  effective  position  angle  of  the  switch's  symmetry  plane. 
No  change  could  be  detected  between  the  position  angle  with  or  without 
the  probe. 

Eecause  the  only  effect  of  the  probe  was  to  change  the  DC  level 
of  the  first  two  months  data,  this  effect  was  subtracted  out  by 
adjusting  the  level  of  these  data  to  agree  w;J_th  the  average  level  of 
the  rest  of  the  year.  Failure  to  do  this  would  have  left  a  step 
function  in  the  data  and  produced  possible  harmonic  contamination  from 
the  spectrum  of  frequencies  contained  in  a  step.   The  only  other 
alternative  would  have  been  to  discard  the  data  and  no  other  reason 
to  do  this  could  be  found. 

It  seems  fairly  certain  that  a  good  portion  of  the  DC  signal  is 
due  to  ground  radiation.   From  equation  2.l6  with  (3  equal  to  zero  we 
obtain  the  expression  for  the  Q  Stokes  parameter,  given  an  axially 
symmetric  universe.   This  parameter  has  a  DC  component  given  by 

Q™  =  (T  -  T  )    cos  6  (1  -  |  sin  $)   ,  (k.l) 

DC   x  w   a'max  2 

where  6  is  the  declination  and  $  is  the  angle  between  the  rotation 

axis  of  the  earth  and  the  symmetry  axis  of  the  universe.  For  (3  equal 

to  ^5  degrees  we  get  the  U  Stokes  parameter  for  which  the  predicted 

DC  value  is  zero.  The  experimentally  measured  values  for  these  two 

quantities  are 

QDC  =  -  .67  ±  .Ik   mdeg  K, 


U__  =  -  .83  ±    .1'+  mdeg  K. 


(h.2) 
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It's  interesting  to  note  that  the  U  component  which  is  supposed  to  be 
zero  is  larger  than  Q,,  the  non-zero  component  predicted  by  the  model. 
Taking  this  with  the  fact  that  the  resulting  vector  points  right  at 
the  elevator  tower  argues  for  this  being  a  local  effect  due  to  ground 
radiation.  Yet  the  possibility  of  a  nonzero  contribution  from  the 
galaxy  or  an  asymmetrically  expanding  universe  cannot  be  ruled  out. 
The  estimates  of  the  ground  radiation  are  nowhere  near  accurate  enough 
to  allow  separating  out  the  ground  contribution  in  an  effect  this 
small.   For  this  reason  the  result  for  Q,  in (k.$  must  be  taken  as  an 
upper  limit  on  the  DC  polarization  from  an  asymmetrically  expanding 
universe: 

QU     <   .85  mdeg  K,with  90$  confidence. 

This  ends  the  analysis  of  the  DC  part  of  the  signal.  Next  I 
will  subtract  the  respective  DC  component  from  each  Q,  and  U  value 
and  investigate  the  result  for  evidence  of  a  cosmological  polarization. 
The  result  of  the  subtraction  is  given  in  Table  k.l. 

Wow  I  must  ask  whether  or  not  there  is  any  significance  to  what 
is  left  after  the  DC  contribution  is  subtracted?  This  poses  a  statistical 
question.  To  be  answered,  the  statistics  of  the  data  must  be  investigated. 
The  Q,  and  U  parameters  in  Table  k.l   are  averages  of  samples  drawn  from 
normally  distributed  populations  whose  means  are  the  true  signal  in 
each  bin.   Because  the  samples  are  large  (on  the  order  of  several 
hundred  members),  the  population  variance  can  be  replaced  by  the  sample 
variance  defined  as 


TABLE   4.1 
STOKES  PARAMETERS  AFTER  SUBTRACTION  OF  DC   POLARIZATION 
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Q 


T   -  T 
NS   EW 


\ 


U 


Experimental  error  in  Q, 


TNE,SW~TNW,SE 


o       =  Experimental  error  in  U 


N. 


Q 


N 


U 


=  Number  of  measurements  of  Q,  in  average 


=  Number  of  measurements  of  U  in  average 


Bin 


Q 


Folded  in  Solar  Bins 


Q 


N„ 


U 


u 


N 


u 


1 

O.98 

0.6l 

1034 

-O.19 

0.60 

1035 

2 

-0.40 

0.6l 

1044 

-0.31 

0.6l 

1047 

3 

0.16 

0.60 

1051 

-0.24 

0.60 

1054 

4 

-0.24 

0.57 

1053 

0.62 

0.59 

1055 

5 

1.27 

O.58 

1050 

0.22 

O.58 

1050 

6 

0.16 

0.59 

1055 

-0.51 

0.59 

1048 

7 

-0.54 

O.63 

1010 

0.52 

0.62 

1012 

8 

-0.94 

0.65 

936 

0.54 

O.65 

928 

9 

1.32 

0.72 

635 

-0.71 

0.67 

636 

10 

-0.33 

O.87 

445 

0.45 

■  O.85 

447 

11 

-0.U6 

1.10 

283 

-1.84 

1.10 

279 

12 

-0.57 

1.27 

224 

-1.07 

I.27 

217 

13 

-0.1+6 

1.24 

211 

-0.13 

1.24 

214 

l4 

-0.53 

1.10 

274 

-0.73 

1.11 

274 

15 

-1.82 

0.93 

437 

-O.63 

0.95 

436 

16 

-0.68 

0.80 

578 

0.01 

0.79 

568 

17 

-0.5k 

O.69 

830 

O.89 

0.70 

833 

18 

0.82 

0.71 

936 

0.10 

O.65 

933 

19 

-0.73 

O.67 

933 

0.25 

O.67 

918 

20 

0.10 

0.64 

9S3 

0.57 

0.65 

962 

21 

O.9O 

0.61 

975 

-1.34 

0.61 

975 

22 

0.16 

0.60 

1000 

0.4D 

0.60 

1003 

23 

0.72 

0.60 

1037 

0.33 

0.59 

1034 

24 

-1.25 

0.62 

1038 

0.16 

0.60 

1039 

TABLE  4.1  (continued) 
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Bin      Q        Q 
Folded  in  Sidereal  Bins 


N 


Q 


U 


Ju 


Nu 


1 

0.0 

O.67 

779 

O.67 

0.66 

776 

2 

-0.24 

O.67 

813 

0.24 

0.65 

805 

3 

0.44 

0.72 

813 

1.53 

0.68 

811 

k 

-1.62 

0.66 

809 

1.13 

0.67 

813 

5 

0.65 

0.65 

810 

0.47 

0.67 

815 

6 

-0.13 

0.64 

809 

a.  98 

0.67 

800 

7 

-0.15 

0.66 

769 

0.21 

0.66 

773 

8 

-i.4o 

0.70 

730 

0.06 

0.67 

718 

9 

1.14 

0.72 

706 

-0.17 

0.74 

704 

10 

-0.09 

0.74 

681 

0.03 

0.74 

681 

li 

0.83 

0.74 

664 

-0.26 

0.72 

663 

12 

0.09 

0.73 

648 

-0.29 

0.74 

64l 

13 

1.06 

0.80 

646 

o.o4 

0.80 

64l 

14 

0.75 

0.82 

643 

0.72 

0.80 

642 

15 

I.89 

0.78 

708 

-0.64 

0.77 

710 

16 

0.19 

0.68 

838 

-0.95 

0.68 

847 

17 

0.14 

0.62 

1003 

-0.56 

0.62 

989 

18 

-0.23 

0.62 

1024 

-0.08 

0.62 

1014 

19 

-O.87 

0.60 

1022 

0.18 

O.58 

1021 

20 

0.15 

0.63 

1008 

-0.69 

O.63 

1027 

21 

1.10 

0.69 

842 

-0.01 

0.71 

839 

22 

-0.55 

0.72 

762 

-1.66 

0.72 

773 

23 

-0.79 

0.73 

732 

-1.45 

0.73 

734 

24 

-1.03 

0.70 

773 

o.i4 

0.70 

760 

Gk 


N       ^2 


°2  =  A.^\-v 


1=1  (h.k) 

=  N  ill  (Xi  "  X>     ' 

It  should  be  noted  that  this  is  not  the  variance  of  the  original 
measurements  and  is  therefore  not  a  fair  estimate  of  the  experimental 
error.  When  I  performed  the  subtraction  between  points  90  degrees 
apart  in  position  in  order  to  subtract  out  the  switch  offset,  I  formed 
a  new  data  set  by  averaging  adjacent  points  in  the  old.  When  each 
point  is  assigned  its  characteristic  sign  dependent  upon  polarimeter 
position,  the  process  may  be  written  as 

Pi'  -  Pi+  Pi  ±  1,  C*.5) 

2 

where  P  is  a  given  Stokes  parameter.   By  the  propagation  of  errors  one 

can  show  that  the  variance  of  the  new  data  set  is  related  to  the  old 

by 

2    GP  2 
Pi        2 

Thus  in  terms  of  the  original  measurements,  Equation  h.k   is  an  estimate 

of  the  variance  divided  by  two;  therefore,  I  must  multiply  by  two 

to  get  back  a  true  estimate  of  the  scatter  in  the  data. 

There  still  remains  one  step  left  in  order  to  determine  the  error 

in  the  measurements  as  they  are  binned.  As  the  data  is  folded,  five 

consecutive  readings  from  each  day's  data  are  averaged  into  a  given  bin. 

Because  of  the  method  of  subtraction,  these  five  points  are  derived 

from  six  independent  points  in  the  original  set  of  measurements. 

Because  a  whole  day  intervenes  between  each  successive  group  of  five 

measurements  in  a  particular  bin,  these  sets  of  five  are  independent 
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and  averaging  over  all  the  points  in  a  bin  is  just  equivalent  to  averaging 

together  the  averages  of  each  day's  contribution.   This  daily  average 

can  be  written  a.s 

P.+   2  P.    .+   2  P.    0+   2  P.    0+   2  P.    ,+   P.    _ 
p.        _     i  i+l  i+2  i+3  1+4       i+5  ,k  7n 

FAvg  "  10  ^'f} 

and  has  a  variance  given  by 

2 
*p.     =   ^i     ,  '  (4.8) 

Avg    556 

2 
where  ap   is  the  variance  of  the  original  measurements.   Just  computing 

i 
Equation  4.4  for  the  bin  in  question  and  dividing  by  the  number  of 

points  would  give  _ 

oS,   -  \  •  (4.9) 

Avg    5 

Therefore,  in  order  to  compute  the  variance  of  the  mean  for  each  of  my 

o 

bins,  I  must  compute  a from  Equation  4.4,  multiply  by  2/N  to  get 

2  i 

o"   ,   and  finally  multiply  by  5.0/5-56  to  correct  for  the  binning 

P. 

procedure. 

In  fitting  a  DC  level  to  each  set  of  Q,  and  U,  I  used  the  method  of 
least  squares,  which  is  to  say,  I  minimized  the  function 

24  (P.-  P)2 

F(P)  =  2   -i§ ,  (4.10) 

i=l   a~ 

i 

2 
where  o_  refers  to  the  corrected  variance  of  the  mean  described  above. 

P. 

l 

If  P  is  the  true  mean  of  the  population  in  each  bin,  then  F(P)  is 

p 
distributed  as  x  with  24  -  1  =  23  degrees  of  freedom.  The  hypothesis 

that  P  is  the  mean  of  each  bin  is  just  the  statement  that  the  polariza- 
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tion  is  isotropic  for  that  particular  Stokes  parameter.  Comparing  the 

actual  minimized  value  of  F(p)  to  the  x  distribution  function  for 

23  degrees, tests  whether  the  isotropic  hypothesis  is  valid.   The 

values  of  F(P)  .   for  each  of  the  fits  to  the  Q  and  U  Stokes  para- 
v  'mm 

meters  and  the  probabilities  that  each  value  or  greater  would  be 
given  by  a  set  of  random  numbers  are: 


! 


SOLAR  BINS 


U: 


2 
X 
30.60 

17.93 


Probability 
15$ 

1% 


SIDEREAL  BINS 

2 
X 
32.98 

27.07 


Probability 
10$ 


Except  for  the  U  parameter  binned  in  solar  time,  the  probabilities  of 
the  parameters  having  a  completely  random  character  is  small.  This 
means  that  statistically  the  measured  polarization  is  not  isotropic, 
or  at  least  the  probability  of  its  being  isotropic  on  the  basis  of 
this  data  is  small.   I  must  now  determine  whether  this  anisotropy 
has  any   real  meaning  for  cosmology  or  whether  it  can  be  explained  in 
terms  of  more  local  effects.   In  our  model  of  an  axlally  symmetric 
universe,  in  addition  to  the  DC  polarization,  we  are  interested  in  the 
sinusoidal  variations  of  the  Q  and  U  parameters  which  have  12  and  2k 
hour  sidereal  periods.  Because  I  am  dealing  with  a  non-uniform,  finite 
data  string,  I  must  also  be  concerned  with  contamination  of  the 
sidereal  contributions  by  strong  signals  in  the  corresponding  solar 
frequencies. 

The  last  effect  can  arise  in  two  ways.   First,  because  we  are 
dealing  with  finite  data  sets,  adjacent  bins  in  the  power  spectrum 
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U8 
will  not  be  completely  uncorrelated.    A  strong  peak  in  one  bin  will 

enhance  an  adjacent  bin  by  as  much  as  15  percent  of  the  value  of  the 

strong  peak.  This  occurs  because  the  elementary  frequency  bins  are 

2    2 
no  longer  rectangular  in  shape,  but  fall  off  as  sin  x/x  where  x  is 

in  units  of  one  over  the  length  of  the  data  string  in  seconds .  When 

the  data  sets  are  not  uniform  in  variance  over  their  length  (as  when 

data  are  thrown  away),  the  sine  and  cosine  functions  used  to  extract 

the  Fourier  components  from  which  the  power  spectrum  is  derived  are  no 

longer  orthogonal.   This  increases  the  correlation  between  adjacent 

bins.   For  this  reason  we  must  be  careful  that  all  large  solar  effects 

are  screened  out.  This  is  particularly  true  of  the  2k   hour  solar 

and  sidereal  components,  which  are  adjacent  elementary  bins  in  the 

power  spectrum  after  one  year  of  data. 

Another  way  in  which  solar  effects  can  hurt  is  through  seasonal 
effects  modulating  the  daily  solar  cycles.   The  product  of  a  2k   month 
seasonal  variation  with  a  solar  period  gives  a  beat  frequency  equal 
to  the  corresponding  sidereal  period.   To  look  for  these  effects, 
the  best  starting  point  is  to  examine  the  data  folded  into  solar 
bins  for  a  signal  which  would  be  large  enough  to  contaminate  the 
sidereal  results. 

Table  k.2   gives  the  results  of  fitting  the  eight  lowest 
frequency  sine  waves  to  the  data  in  both  solar  and  sidereal  time. 

The  procedure  was  to  minimize  the  function 

p 
2k     (P.-  A^Sinut.-  B„coscot.  ) 


ill 

o. 


P. 

1 
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with  respect  to  the  parameters  A_  and  B„.   From  A  and  B„  one  can 
compute  an  amplitude  and  a  phase  for  the  best  fit  sine  function 


2   „  2XX 


Amplitude  =  (A   +  B  )2  , 


B, 


Phase  =  ARCTAN  (  f ) , 


(4.12) 


Af 


for  each  frequency  component.  Starting  with  the  2k   hour  component,  each 

wave  was  computed,  subtracted  from  the  data,  and  the  residuals  used  to 

2  £ 

compute  ax.  At  each  step  I  compared  the  x  "to  a  random  distribution 

with  two  less  degrees  of  freedom.  The  next  fit  then  uses  the  residuals 

of  the  preceding  step.  Table  k.2   also  gives  the  expected  value  for  each 

frequency  component.   Because  of  the  non-zero  variance  of  each  data  point, 

there  is  an  expected  value  of  the  amplitude  of  each  frequency  component 

given  by  g 

/  2k  sin  cut      r  ,2k  cos   cut.  \2        ,    2k  coscut  sincut       gn 

<Amplitude>    =  i^i __i  UJ    i  )      '   ^±ti £ )    -I 

o±-  o.^  „* 


2 


r2k     coscut.      2,k     sincut.  ,2k     coscut.  sincut.    N2-,2 

CLl  — — ^    tl  — ^  -  (tl  -** i  )   3 

a. 2  a.  «=  ct.2 

11  l 


24     cos  cut.  2k  sincut.      _  2k     coscut.  sincut.      _ 

ill     "^  [(A  -^ri  )2   -    (f=l      ,     0^2         X    )21 

+ 1 1 1 

r2k     coscut.      2k     sincut.         ,2k     coscut.  sincut.    s.2n2 
L2  ,  i     2  ..  i-(S.  i  l        J 

i=l  p —  1=1  P —    1=1  5 

a.  a.  ^  o.  <=■ 

11  l 

(4.13) 

To  test  the  statistical  significance  of  an  amplitude,  one  takes  the 

difference  between  the  amplitude  and  its  expected  value  for  random  data 
and  compares  this  with  the  error  in  the  amplitude. 

Investigation  of  the  solar  amplitudes  shows  only  one  with  a 
significant  value,  the  k  hour  period  in  the  Q  parameter.  This 
frequency  is  well  separated  from  its  sidereal  equivalent  in  the  power 

spectrum,  as  one  needs  only  two  months  of  data  to  resolve  the  two. 


TABLE  4.2 
HARMONIC  ANALYSIS  OF  DATA 
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P  =  Period  in  hours  of  solar  or  sidereal  time. 

T  =  Amplitude  of  harmonic  component  in  mdeg  K. 

a  -  Error  in  amplitude. 

§  =  Phase  of  component  in  degrees  from  Bin  0  (Bin  24). 

a&  =  Error  in  phase. 

v  =  Degrees  of  freedom  in  fit . 

2 

X  =  Chi-square  for  the  fit. 

<T>  =  Expected  value  of  T  based  on  experimental  errors  and  random 
numbers. 


T-<T> 


=     Number  of  sigrna  of  amplitude  above  expected  value 


<T> 


T-<T> 


X 

* 

CTT 

Q  in  Solar 

Bins 

2k       0.33 

0.21 

97.36 

34.35 

21 

26.94 

0.29 

0.19 

12   0.19 

0.20 

-123.61 

62.45 

19 

25.59 

0.29 

-0.50 

8   0.28 

0.20 

-42.31 

42.03 

17 

24.21 

0.29 

-0.05 

6   0.22 

0.20 

-100 . 89 

53.08 

15 

23.25 

0.29 

-0.35 

4.8  0.09 

0.20 

182.42 

129.61 

13 

23.22  - 

0.29 

-1.00 

k       0.60 

0.20 

73.14 

19.40 

11 

14.95 

0.29 

1.55 

3.4  0.23 

0.20 

30.88 

51.09 

9 

14.96 

0.29 

-0.30 

3   0.15 

0.20 

-89.83 

78.45 

7 

14.64 

0.29 

-0.70 

U  in  Solar 

Bins 

2k       0.10 

0.21 

123.20 

109.79 

21 

16.99 

0.29 

-0.90 

12   0.23 

0.20 

-40.31 

49.17 

19 

15.35 

0.29 

-0.30 

8   0.l4 

0.20 

139.9^ 

84.06 

17 

14.45 

0.29 

-0.75 

6   0.19 

0.20 

-163.76 

60.61 

15 

13.92 

0.29 

-0.50 

k.8  0.31 

0.20 

-102.86 

37.02 

13 

11.79 

0.29 

0.10 

k       0.24 

0.20 

-133-50 

48.93 

11 

10.96 

0.29 

-0.25 

3.k   0.08 

0.20 

-124.28 

143.00 

9 

10.92 

0.29 

-1.05 

3   0.47 

0.20 

52.36 

24.81 

7 

5.61 

0.29 

0.90 
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TABLE  4,2    (Continued) 


p 

T 

Om 

§ 

a. 

V 

T 

§ 

Q  in 

Sidereal 

Bins 

24 

0.52 

0.21 

-80.52 

21.79 

21 

12 

0.20 

0.20 

89.61 

58.59 

19 

8 

0.19 

0.20 

-83.75 

61.66 

17 

6 

0.31 

0.20 

-16.64 

37.81 

15 

4.8 

o.4i 

0.20 

8k.  80 

28.58 

13 

4 

0.36 

0.20 

62. 6k 

31.95 

11 

3.4 

0.10 

0.20 

-59.81 

116.81 

9 

3 

0.39 

0.20 

-105.69 

30.13 

7 

U  in 

Sidereal 

Bins 

2k 

0.58 

0.20 

39.24 

20.53 

21 

12 

0.45 

0.21 

22.62 

25.30 

19 

8 

0.33 

0.20 

-1.32 

35.01 

17 

6 

0.30 

0.20 

79.77 

38.67 

15 

k.Q 

0.17 

0.20 

187.52 

67.42 

13 

k 

0.08 

0.20 

40.47 

136.89 

11 

3.k 

O.29 

0.20 

133.51 

40.62 

9 

3 

0.25 

0.20 

-163.63 

46.74 

7 

2 

X 

<r> 

T-<T> 

25.65 

0.29 

1.10 

24.35 

0.29 

-0.45 

23.10 

0.29 

-0.50 

20.89 

0.29 

0.10 

17.11 

0.29 

0.60 

13.68 

0.29 

0.35 

13.63 

0.29 

-0.95 

10.42 

0.29 

0.50 

17.30 

0.29 

1.45 

12.78 

0.29 

0.76 

9.63 

0.29 

0,20 

7.14 

0.29 

0.05 

6.46 

0.29 

-0.60 

6.28 

0.29 

-1.05 

3.85 

0.29 

0.00 

2.61 

0.29 

-0.20 
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The  absence  of  significant  solar  values  with  frequencies  lower  than  this 
indicates  that  contamination  of  sidereal  frequencies  with  solar  signal 

won't  be  a  problem.  The  solar  signal  is  not  zero  as  evidenced  by  the 

2 
high  initial  x  hut  I  would  be  a  little  surprised  if  it  were,  since 

most  environmental  variations  which  can  affect  the  polarimeter  are 

solar  in  nature.   It  appears  that  by  shielding  the  polarimeter  from 

ground  radiation  and  by  throwing  away  data  when  the  sun  was  above  the 

ground  shield,  I  was  able  to  avoid  a  large  amount  of  contamination. 

Before  discussing  the  data  further,  there  is  one  feature  of  a  fit 

of  this  type  that  should  be  noted.   Ordinarily  when  we  fit  sines  and 

cosines  to  a  set  of  numbers  with  uniform  variance,  the  sines  and  cosines 

are  orthogonal  functions  over  the  set.   In  my  case  the  variances  of 

the  points  are  not  equal  but  depend  on  the  number  of  measurements  per 

point  and  the  actual  scatter  in  those  measurements.   To  insure  that  I 

wasn't  adding  any  systematics  to  the  data  in  the  fitting  process,  I 

compared  the  results  obtained  above  to  the  results  one  would  obtain 

assuming  a  constant  variance  across  the  data  set.  For  each  set  I  chose 

as  a  variance  the  average  variance  per  point  of  the .corresponding 

Stokes  parameter  folded  into  sidereal  bins.  The  results  are  given  in 

Table  k.3.      The  most  notable  differences  are  in  the  solar  frequencies. 

Here  the  points  closest  to  noon,  which  were  deweighted  because  of 

fewer  measurements  (due  to  solar  contamination),  were  given  equal 

weight.   It  is  also  these  points  which  are  most  affected  by  the  sun 

and  therefore,  it  is  not  surprising  that  the  solar  amplitudes  become 

somewhat  larger.   The  effect  on  the  sidereal  data  is  small.   In  both 

cases  the  Chi -squares  seem  to  follow  the  same  general  progression  as 
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TABLE  4.3 
HARMONIC  ANALYSIS  OF  DATA  ASSUMING  UNIFORM  VARIANCE 
P  T  $  v 

Q  in  Solar  Bins 


24 
12 

8 

6 

4.8 

4 

3.4 

3 


0.43 

93.32 

21 

27.44 

0.31 

-114.01 

19 

24.97 

0.22 

-13.12 

17 

23.68 

0.17 

-109.24 

15 

22.92 

0.12 

-117.86  t 

13 

22.58 

0.59 

75.56 

11 

13.73 

0.14 

9.71 

9 

13.20 

0.12 

-93.50 

7 

12.86 

U  in  Solar  Bins 

24 
12 

8 

6 

4.8 

4 

3.4 

3 

Q,  in  Sidereal  Bins 

24 
12 

8 

6 

4.8 

4 

3.4 

3 

U  in  Sidereal  Bins 

24 
12 

8 

6 

4.8 

4 

3.4 

3 


0.29 

121.77 

21 

21.17 

0.39 

-46.24 

19 

17.26 

0.24 

140.08 

17 

15.72 

0.08 

-164.68 

15 

15.55 

0.33 

-105.00 

13 

12.71 

0.12 

181.05 

11 

12.31 

O.21 

-93.80 

9 

11.15 

0.45 

66.11 

7 

5.96 

O.56 

-83.84 

21 

27.38 

O.25 

82.76 

19 

25.78 

0.22 

-83.32 

17 

24.58 

0.34 

-8.13 

15 

21.66 

0.39 

92.11 

13 

17.86 

0.35 

56.00 

11 

14.78 

0.08 

-50.59 

9 

14.63 

O.38 

-111.56 

7 

10.85 

0.59 

34.81 

21 

19.29 

0.46 

20.38 

19 

13.68 

0.35 

0.05 

17 

10.55 

0.34 

78.64 

15 

7.63 

O.18 

184.40 

13 

6.81 

0.10 

61.80 

11 

6.56 

0.34 

189.77 

9 

3.63 

0.21 

-157.64 

7 

2.46 
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in  the  weighted  analysis,  indicating  that  no  systematics  were  introduced 
due  to  the  non-uniform  weighting. 

Turning  to  the  sidereal  frequencies,  there  are  three  significant 
amplitudes,  both  2k   hour  components  and  the  12  hour  component  for  the 
U  parameter.   The  largest  of  these  is  the  2k   hour  component  of  U.  The 
positive  maximum  of  the  fitted  sine  wave  falls  at  a  right  ascension  of 
2.88  ±  1.37  hours  which  is  about  two  hours  a$*ay  from  the  plane  of  the 
galaxy  in  the  anticenter  region.  The  negative  minimum  falls  in  the 
region  of  the  North  Galactic  spur  at  1^.88  ±  1.37  hours.  Looking  at 
the  surveys  of  galactic  continuum  polarization  at  lower  frequencies, 
these  regions  seem  to  exhibit  significant  amounts  of  polarization. 
Looking  also  at  the  plot  of  the  U  data  in  Figure  k.k,   we  note  that  in 
addition  to  the  high  values  around  bins  3  and  k   arid  the  low  values 
around  bins  15  and  l6,  there  is  another  low  region  extending  from  bin 
20  to  bin  23,  or  19.5  to  22.5  hours  in  right  ascension.  Again  looking 
at  a  galactic  map,  I  find  that  the  Cygnus  arm  of  the  galaxy,  the 
hottest  region  to  pass  through  my  beam,  falls  at  20.5  hours.   Because 
I  cut  the  plane  of  the  galaxy  in  two  places  roughly. 12  hours  apart 
in  right  ascension,  it  isn't  too  difficult  to  imagine  a  12  or  2k   hour 
component  being  generated  if  any  polarized  component  existed  in  these 
regions. 

From  a  galactic  polarization  survey  at  21  cm   I  found  that 
individual  measurements  of  100  to  200  millidegrees  were  obtained  in  the 
regions  I  cover.  Applying  a  spectral  index  of  a  -2.7  to  this  intensity, 
in  order  to  extrajjolate  the  reading  to  3.2  cm,  results  in  a  possible 
polarization  temperature  of  from  .6  to  1.2  millidegrees  Kelvin. 


7k 


Although  this  arm  waving  proves  little,  it  did  convince  me  that  I 
had  better  look  into  polarized  galactic  radiation  in  detail. 

k.k     SUBTRACTION  OF  THE  GALAXY  COIWEIBUTIOW 

In  any  type  of  extrapolation,  one  would  like  to  use  data  which  has 
parameters  as  close  to  those  of  the  data  being  corrected  as  possible. 
For  a  radio  survey  of  the  galaxy,  the  two  most  important  parameters 
are  frequency  and  angular  resolution.   Of  these,  frequency  has  the 
edge  in  importance  when  discussing  a  polarization  survey. 

There  are  two  ways  in  which  frequency  enters  into  extrapolating 
polarization  measurements.  Most  directly  it  must  be  considered  in 
terms  of  the  spectrum  of  the  radiation  one  is  measuring.  Most  of  the 
galactic  polarized  emission  is  thought  to  be  non-thermal  synchrotron 
radiation  with  a  spectral  index  of  -2.5  to  -2.9.       By  spectral 
index  we  mean  the  quantity  (3  in  the  expression 

1     1 
for  scaling  the  temperature  of  radiation  at  different  frequencies.   If 
the  radiation  follows  this  simple  relation,  knowing  the  spectral  index 
determines  the  magnitude  of  the  polarization  at  one  frequency,  given 
the  magnitude  at  another . 

The  other  aspect  of  polarized  emission  in  the  galaxy,  the  angle, 
also  depends  on  frequency  through  the  mechanism  of  Faraday  rotation. 
The  electron  concentration  in  the  galaxy  together  with  the  galactic 
magnetic  field  act  to  produce  rotation  of  linear  polarization  in  much 
the  same  way  a  magnetic  field  acting  on  a  ferrite  produces  rotation 
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in  the  Faraday  rotation  switch.  This  effect  is  frequency  dependent, 

53 
satisfying  a  relation  of  the  form 


2 
*  =  ill  +  R  a   , 

T    o    m    ' 

where  R  is  called  the  rotation  measure  and  is'  given  in  units  of 
m 


(^•15) 


radians  per  meter  squared. 


and  •{/     are  the  angles  of  the  radiation 


before  and  after  rotation.   In  terms  of  the  number  of  electrons  and  the 

longitudinal  magnetic  field  along  the  line  of  sight,  the  rotation 

53 
measure  may  be  written  as 


R     =  8.1  x  10 5    P  N  BTdL     RADIANS/]/ 
m  -J  L  ' 


(h.l6) 


Because  N  and  B  are  not  constant  with  direction  but  vary  in 

h 


a  complex  way,  one  cannot  use  them  as  fitting  parameters  to  align  two 
surveys  at  different  frequencies.   One  needs  at  least  two  surveys  at 
different  frequencies  to  make  a  correction  for  a  third.  An  even  better 
situation  is  to  have  three  surveys  to  fit  the  two  parameters  in  Equation 
U.10.   This  would  give  an  extra  degree  of  freedom  with  which  to  test 
independently  for  a  good  fit.  There  are  three  surveys  in  the  literature 
at  Uo8,  6l0  and  1^07  MHz  with  which  to  attempt  to  determine  the  rotation 


measure  in  my  region  of  interest 


^9,5^,55 


Although  the  resolution  of  all  three  surveys  is  approximately  the 
same,  2  of  arc,  they  don't  cover  the  same  points.   Because  of  this 
and  because  I  have  to  wind  up  with  a  result  which  is  averaged  over  my 
beamwidth,  I  decided  to  convolve  the  surveys  with  my  antenna  pattern 
before  attempting  to  fit  Equation  U.10.   By  doing  so  I  felt  that  any 
systematics  in  individual  readings  or  sparse  coverage  in  a  region  for 
a  particular  survey  would  be  averaged  out  in  the  large.  Table  k.h 
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TABLE  4.4 

BINNED  SUEVEYS  AT  LOWER  FREQUENCY 

BINNED  408  MHz 

SURVEY 

Bin 

Q 

U 

T 

CTT 

Angle 

AAngle 

1 

0.2551 

-0.8107 

0 . 8499 

0.0206 

143.7336 

O.6947 

2 

1.21+17 

-0.0071 

1.2417 

0.0215 

179.8353 

0.4954 

3 

0.1+868 

0.2121 

0.5310 

0.0210 

11.7738 

1.1343 

k 

0.3362 

-0.1617 

0.3731 

0.017O 

167.1560 

1.3049 

5 

-O.41+97 

-0.2569 

0.5179 

0.0145 

104.8694 

0.8029 

6 

-O.7188 

-0.4034 

0.8242 

o.oi48 

104.6497 

0.5137 

7 

-O.7658 

-0.2550 

0 . 8072 

0.0173 

99 . 2100 

0.6157 

8 

-0.021+0 

-0.5514 

0.5520 

0.0223 

133.7520 

1.1585 

9 

0.1274 

-0.6139 

0.6270 

0.04*41 

140.8636 

2.0131 

10 

O.3056 

-0.5857 

0.6606 

0.0330 

148.7786 

1.4292 

11 

0.5183 

-0.4465 

0.6841 

0.0221 

159.6279 

0.9245 

12 

0.1812 

-0.2787 

0.3325 

0.0196 

151.5155 

1.6904 

13 

-0 . 1+166 

-0.0449 

0 . 4190 

0.0312 

93.0771 

2.1354 

III 

0.0855 

0.6023 

0.6083 

0.0159 

40.9592 

0.7475 

15 

-0.1455 

0.5694 

0.5877 

0.0115 

52.1664 

0.5588 

16 

-O.1987 

0.5723 

0.6058 

0.0121 

54.5732 

0.5702 

17 

0.2000 

0.6266 

0.6577 

0.0131 

36.1509 

0.5706 

18 

0.3639 

0.5715 

0.6775 

0.0148 

28.7552 

0.6268 

19 

0.2257 

0.5482 

0.5928 

0.0160 

33.8108 

0.7746 

20 

-O.I929 

0.3432 

0.3937 

0.0151 

59.6691 

1.0996 

21 

-0.2014 

0.1726 

O.2652 

0.0176 

69.6953 

I.8996 

22 

0.1154 

0.0236 

0.1178 

0.0230 

5.7685 

5.5845 

23 

0.0574 

-0.0345 

O.0670 

O.0167 

164.4988 

7.1199 

24 

-0.3651+ 

0.1052 

O.3802 

0.0186 

81.9654 

1 . 4010 

BINNED  610  MHz  SURVEY 


Bin 

1 
2 

3 

4 

5 
6 
7 
8 

9 
10 
11 
12 

13 
11+ 
15 
16 
17 
18 


Q 

0.1924 
O.2718 
0.2162 
0.1281 
-0.0933 
-0.2577 
-O.2955 
-O.1610 
-O.0889 
-0.7657 
-0.3553 
-O.1789 
-0.0441 
-0.04L3 
-0.0885 
-0.1267 
-0.5166 
-o.6i84 


u 

-0.0906 

-0.0236 

0.0383 

0.1774 

0.1882 

-0.0665 

-0.0299 

-0.1064 

-0.1318 

-0.1950 

-0.2296 

-0.3978 

-0.1889 

0.1135 

0.1691 

0.1612 

0.1374 

0.3182 


0.2127 
0.2728 
0.2196 
0.2189 
0.2101 
0.2661 
0.2971 
0.1930 
0.1590 
0.7901 
0.4231 
0.1+362 
0.1940 
0.1208 
0.1908 
0.2051 
0.5345 
0.6954 


T 

0.006l 
0.0067 
0.0057 
0.0047 

0.0039 
0.0039 

0.0043 
0.0050 
0.0057 
0.0057 
0.0043 
0.0034 
0.0032 
0.0031 
0.0032 
0.0034 
0.0037 
o.oo4i 


Angle 

167.3881 

177.5216 

5.0284 

27.0834 

58.1868 

97.2307 

92 . 8929 

106.7299 

117.9949 

97.1432 

106.4366 

122.8918 

128.4281 

54.9980 

58.8083 

64.0847 

82.5539 

76.3853 


AAngle 

0.8248 
0.7058 
0.7452 
o.6i44 

0.5368 
0.4185 

0 . 4i4i 
0.7448 
1.0280 
0.2085 
0.2926 
0.2224 
0.4664 
0.7461 
0.4837 
0.4763 
0.1973 
0.1702 


TABLE  4.4  (continued) 
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BINNED  6lO  MHz  SURVEY  (continued) 
Bin      Q,         U         T 


Angle    AAngle 


19 

-0.7566 

0.6826 

I.0190 

0.0046 

68.9721 

0.1293 

20 

-0.4072 

0.4526 

0.6088 

0.0043 

65.9895 

0,2016 

21 

-0.1798 

0.2445 

0.3035 

o.oo4i 

63.1695 

0.3915 

22 

0.0294 

0 . 5464 

0.5472 

0.0042 

43.4574 

0.2201 

23 

o.i44i 

0.3553 

0.3835 

O.OQ46 

33.9612 

0.3403 

24 

0.0594 

O.O193 

O.0625 

0.0054 

8.9797 

2.4979 

BINNED  1407  MHz  SURVEY 


Bin 


Q 


U 


JT 


Angle    AAngle 


1 

0.0524 

-0.0280 

0.0594 

0.0058 

165.9682 

2.8037 

2 

0.1347 

0.0405 

0.1407 

O.OO58 

8.3572 

1.1818 

3 

o.ii4o 

0.0902 

0.1454 

O.OO58 

19.1729 

1.1340 

4 

0.0185 

0.0619 

0.0646 

0.0055 

36.6985 

2.4473 

5 

-0.0260 

0.1522 

0.1544 

O.OO56 

49.8548 

1.0467 

6 

-0.0926 

0.1388 

O.1668 

0.0055 

61.8627 

0.9396 

7 

-0.1452 

0.0527 

0.1544 

0.0053 

80.0338 

0.9861 

8 

-0.1256 

0.0104 

O.1260 

0.0055 

87.6273 

1.2446 

9 

-0.0540 

-0.0008 

0.0540 

0.0054 

90.4255 

2.8782 

10 

-O.OO9O 

0.0070 

0.0114 

0.0053 

70.9943 

13.3585 

11 

0.0195 

-0.0063 

0.0205 

0.0053 

170.9830 

7.3773 

12 

0.0054 

-0.0365 

O.0368 

0.0052 

139.2055 

4.0506 

13 

0.0113 

-0.0537 

0.0549 

0.0049 

140.9456 

2.5499 

14 

0.0843 

-0.04-10 

0.0937 

0.0047 

167.0171 

1.4336 

15 

0.1101 

-0.0382 

O.1166 

o.oo44 

170 . 4431 

1.0888 

16 

0.0875 

0.0158 

O.0889 

0.0044 

•5.1105 

1.4255 

17 

-0.0199 

0.0815 

O.0839 

o.oo44 

51.8526 

1.5014 

18 

-O.0685 

0.0232 

0.0724 

o.oo44 

80.6350 

I.7262 

19 

-0.0701 

0.0202 

O.0729 

o.oo46 

81.9476 

1.7981 

20 

-O.0522 

-0.0112 

0.0534 

o.oo48 

96.0603 

2.5859 

21 

-0.0465 

-O.O694 

0.0835 

0.0051 

118.0791 

1.7480 

22 

-0.0303 

-O.O769 

O.0826 

0.0053 

i2i+.2484 

1.8538 

23 

0.0190 

-0.0859 

0.0880 

0.0050 

141.2228 

I.6261 

24 

0.0353 

-O.O92O 

O.O986 

0.0051 

145.4948 

1.4855 
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gives  the  results  of  the  convolution,  where  the  convolved  data  have  been 
binned  in  exactly  the  same  way  as  my  actual  data. 

In  attempting  to  fit  the  angles  of  these  data  to  Equation  U.10,  I 
fixed  one  of  the  three  angles  and  tried  the  fit  with  various  multiples 
of  180  degrees  added  to  the  other  two  angles.  '  In  this  way  I  sought  to 
account  for  possible  rotations  of  greater  than  180  degrees  between 

frequencies.   I  then  selected  the  best  fit  fpr  each  bin  based  on  the 

2  2 

X  and  have  displayed  the  results  in  Table  4.5.   In  most  cases  the  x 

is  huge,  though  it  should  be  less  than  one  for  a  reasonable  fit.  An 

explanation  is  that  the  coverage  in  each  individual  survey  is  not 

good  enough  to  allow  accurate  averaging.  Coverage  is  an  important 

factor  in  that  the  electron  concentrations  and  fields  change  significantly 

on  a  small  scale,  at  times  even  less  than  that  of  the  resolution  of  the 

individual  surveys.   If  not  enough  of  the  contributions  are  added,  the 

average  becomes  inaccurate.  Another  possibility  is  that  localized 

Faraday  depolarization  causes  different  areas  at  different  frequencies 

to  be  important  in  the  average  of  a  particular  bin. 

Because  of  the  poor  quality  of  the  fits  obtained  for  the  three 
lower  frequencies,  I  decided  that  the  angles  obtained  from  them  would 
be  useless  for  the  purpose  of  correcting  my  survey  and  that  I  should 
forget  them  entirely  and  concentrate  on  working  with  the  amplitudes. 
I  felt  that  if  there  were  a  reasonable  correlation  between  the  amplitudes 
from  my  survey  and  the  survey  at  the  next  lowest  frequency,  I  could 
fit  the  amplitudes,  using  the  spectral  index  as  a  parameter. 

Amplitudes  obey  Rayleigh  rather  than  Gaussian  statistics.  The 
correlation  coefficient  and  its  associated  expectation  values  for 
Rayleigh  statistics  are  given  by 
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TABLE  k.5 
RESULTS  OF  FITTING  A  ROTATION  MEASURE  TO  LOWER  FREQUENCY  SURVEYS 


Bin               *M  AEM  '  ^ 

1  -32.82  0.003  ik.lG 

2  63.18  0.001  1.92 

3  63.22  0.003  21.kk 
k  -33.73  0.006  17.39 

5  65.51  0.002  69.O7 

6  52.62  0.001  301.79 

7  -31.10  0.001  8.17 

8  1.62  0.003  0.53 

9  -30.06  0.012  1+.31 

10  66.11  0.007  0.03 

11  -38. 94  0.003  3.^7 

12  -50.80  0.009  15.30 

13  -33.48  0.012  3.^9 
Ik  -1+2. 67  0.002  28.85 

15  -U2.3^  0.001  18.28 

16  -U2.63  0.001  0.80 

17  -13.23  0.001  0.08 

18  -3^.20  0.001  150.27 

19  -33.56  0.002  7.48 

20  62.61  0.001+  34.27 

21  -20.73  0.008  1.24 

22  -23.19  0.022  0.02 

23  70. 47  0.019  O.36 

24  35.83  0.005  0.02 
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N 

s  (T.J1-  ^r^)  (t2x-  <r2> 
c^1'2')  =   .4292  n  5"  o~       ' 

<C(l,2)>  =  0,  (H.17) 

2   i 
<C(1,2)   =  N 


For  the  correlation  between  the  3.2  and  21  cm  data,  I  find  that 
C(9370,l407)  =  .1663.  Comparing  this  with  the  expected  error  for  2h 
points  shows  that  this  is  only  about  a  .8  sigma  effect  and  not  a  very 
convincing  correlation. 

In  spite  of  the  low  correlation,  I  still  carried  out  the  fit, 
feeling  that  though  I  could  not  use  it  for  purposes  of  subtraction,  I 
could  at  least  ask  whether  or  not  the  scale  difference  between  the  two 
sets  of  measurements  was  reasonable.  The  value  of  the  spectral  index 
which  minimizes  the  Chi-square  is  -2.^9,  which  doesn't  seem  too  un- 
reasonable considering  the  possibility  of  Faraday  depolarization.   This 
confirms  my  back  of  the  envelope  calculation  and  indicates  that  on  the 
basis  of  the  21  cm  data  I  should  be  seeing  an  effect  at  3.2  cm.   To 
find  the  source  of  the  lack  of  correlation,  I  plotted  the  scaled 
amplitudes  of  the  two  sets  of  data  in  Figure  4.5.  The  biggest 
discrepancy  is  the  lack  of  a  significant  contribution  at  21  cm  between 
19  and  23  hours  of  right  ascension.   In  this  area  is  the  Cygnus  arm  of 
the  galaxy,  the  hottest  region  that  falls  in  the  beam. 

Going  back  to  Table  h.k,   there  doesn't  appear  to  be  a  contribution 
in  this  region  at  408  MHz  either,  yet  looking  at  a  survey  at  the  same 
frequency  and  an  8  bearawidth,  there  is  an  indication  that  an  effect  in 
fact  does  exist.    In  this  region  especially,  where  the  galactic 
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magnetic  fields  are  complex  and  are  directed  along  the  line  of  sight, 
one  would  expect  Faraday  depolarization  and  problems  of  coverage  to  be 
important . 

I  feel  that  there  is  strong  evidence  for  a  contribution  to  the 
3.2  cm  polarization  from  the  galaxy.  Unfortunately,  the  data  available 
at  other  frequencies  does  not  seem  capable  of  supplying  corrections. 
In  all  fairness  I  can  understand  why  they  wouldn't,  given  the  complex 
nature  of  the  emission  regions  and  magnetic  fields  in  the  galaxy,  the 
large  gap  in  wavelength  between  them  and  3.2  cm  ,  and  the  purpose  for 
which  these  surveys  was  made,  which  was  to  provide  high  resolution 
coverage  of  the  most  important  regions  of  polarized  radiation  at  their 
frequencies.   In  their  case  the  most  important  regions  were  close  to 
the  galactic  center  and  therefore  their  coverage  was  concentrated  in 
this  area. 


83 


k.5     FINAL  ANALYSIS 

Being  unable  to  subtract  the  galaxy,  I  have  to  accept  the  harmonic 
components  as  they  are  given  in  Table  k.2   and  investigate  the  upper 
limit  they  infer.   From  here  the  problem  lies  in  the  complicated 
predictions  of  the  model  given  by  equation  2.l6.  Rewriting  this 
equation  twice,  once  with  6=0  and  once  with  6=^5  ,  I  produce  the 
predictions  for  the  Q  and  U  parameters,  at  a  declination  6 

Q  =  (T  -  T  )   [sin2ssin6cos6sin(t-Q'  -rr  )  +  sin  §(l-icos  6)sin2(t-a  +n  )], 
^   v  w   a'max  v   °~p  °  ~TT 

U  =  (T  -  T  )    [sin2§cos6sin(t-a  +rr)  +  sin  §sin6sin2(t-0'  )]  , 
v  w   a'max  v    o  v   oy J  ' 

(^.18) 

where  everything  is  expressed  in  terms  of  sine  functions.   These  may  be 
compared  directly  with  my  data.  The  problem  in  doing  this  lies  in  the 
fact  that  the  phases  of  the  experimentally  measured  components  don't 
agree  well  with  the  model.  This  isn't  surprising,  considering  that  most 
of  the  contribution  is  due  to  galaxy  and  not  an  asymmetrically  expanding 
universe.   For  example,  the  measured  phases  of  the  2k   hour  components 
for  the  Q  and  U  parameters,  when  compared  with  Equation  4.13,  predict 
values  for  a       which  differ  by  150.24  degrees,  well  outside  the  exper- 
imental errors . 

This  suggests  that  by  forcing  a  fit  to  the  model,  the  non-agreement 
of  phases  would  help  toward  a  lower  limit.   I  object  to  this  course  for 
two  reasons.   First  the  error  introduced  by  the  phases  may  be  large 
enough  that  any  lowering  of  the  limit  would  probably  be  offset  by  an 
increase  in  the  uncertainty.  More  importantly,  I  feel  that  because  I 
do  have  statistically  significant  amplitudes,  one  can  always  argue  that 


8k 


even  were  they  to  be  predominately  from  the  galaxy,  they  might  be 

masking  a  significant  cosmological  contribution.  A  fairer  estimate  is 

to  compute  a  90%  confidence  upper  limit  for  each  12  and  2k   hour  harmonic 

of  Q,  and  U  and  forgetting  phases,  compute  a  result  in  each  case  for 

the  absolute  value  of  (T  -  T  )    .   This  procedure  assumes  that  the 

v  w   a  max 

phases  in  the  data  fit  the  model  perfectly  and  uses  the  amplitudes 
to  produce  an  upper  limit.  Another  advantage  of  this  method  is  that 
the  statistics  of  the  amplitudes  are  well  known  and  setting  confidence 
levels  is  relatively  easy.   First  we  form  the  quantity, 

where  T  is  assumed  to  consist  of  a  random  component  plus  a  real  signal, 
and  the  expectation  value  is  based  solely  on  random  errors.  This  is 
just  the  signal  power  divided  by  the  noise  power,  a  parameter  which 

has  a  non-central  Chi- square  distribution  for  two  degrees  of  freedom. 

Gk 
Consulting  the  appropriate  probability  graph  for  this  distribution, 

I  find  the  following  upper  limits 

T  <  .75  mdegK 

T  <  .38  mdegK, 

T  <  .81  mdegK  (4.20) 

T  <  .67  mdegK,  with  90%  confidence. 

Inserting  these  into  the  expressions  for  Q  and  U  along  with  the 
declination  of  the  observation,  ^0.35  degrees,  one  obtains 


2k   hr 

For  Q 

12  hr 

2k   hr 

For  U 

12  hr 
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For  Q:   |(T  -  T  )    I  <  1.6l  mdegK, 

'  W   a  max  (4.21) 

U:   I (T  -  T  )    I  <  1.31  mdegK,  with  90$  confidence. 

We  must  also  consider  the  possibility  of  $  being  very  small.   In 

this  case  the  cosmological  contribution  would  be  buried  in  the  DC 

contribution.   Assuming   $  =  0  and  using  the  upper  limit  on  the  DC 

level,  I  find  from  Equation  4.1 

t 

|(T  -  T  )    I  <  1.47  mdegK,  with  90$  confidence.        (4.22) 
1  v  w   a'max '  3  y  '  v 

Taking  the  largest  of  the  three,  I  find  for  an  upper  limit  on  the  mag- 
nitude of  the  polarization  of  the  microwave  background: 

<  1.6l  mdegK,  with  90$  confidence.  (4.23) 

In  terms  of  a  2.7  degree  microwave  background,  this  amounts  to  .06$ 
polarization.   Inserting  this  result  into  Equations  1.19  we  obtain  for 

e  the  asymmetry  parameter  and  A  h  /h  : 

-4 
|e|  <  8.15  x  10 

<^) 

-r—  <  6.41  x  10"^, 
h   ~  ' 

o 

where  I  have  assumed  the  nominal  values  n  =10  /cm  and  h  =  50km/sec 

o      '         o       ' 

Mpc   .   It  should  be  remembered  that  Equations  1.19  were  derived  assuming 
that  l(t),  the  ionized  fraction,  was  equal  to  one  for  all  times.  In 
general  this  is  an  unreasonable  assumption,  but  Rees  has  shown  that  in 
the  event  of  a  reheating  phase  at  a  redshift  of  z  >  7  (due  perhaps  to 

galaxy  formation),  that  the  solution  applies  and  defines  the  upper  limit 

42 
on  the  expansion  asymmetry. 

In  the  event  that  there  was  no  reheating  between  recombination  at 
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z  «  1000  and  the  present,  the  power  anisotropy  measurements  imply  a 
polarization  smaller  than  that  given  in  4.18,  and  in  this  case 
provide  the  upper  limit. 

k.6     FUTUEE  DIRECTIONS 

As  in  any  enterprise,  there  is  always  room  for  improvement  and 
this  one  is  no  exception.  The  two  most  glaring  limits  on  this  exper- 
iment are  the  galactic  contribution  and  the  effects  of  ground  radiation. 
In  principle, the  latter  could  be  solved  by  improving  the  ground  shield 
and  by  using  an  antenna  with  somewhat  better  sidelobe  characteristics. 
In  the  case  of  the  galactic  contribution, there  doesn't  seem  to  be  any 
reasonably  clean  way  of  correcting  for  it  at  3.2  cm.   For  this  reason, 
I  feel  the  best  solution  is  to  try  a  frequency  where  the  expected 
contribution  from  the  galaxy  is  less  than  the  expected  limit  from 
receiver  noise.   For  example,  at  1.6  cm  and  assuming  a  spectral  index 
of  -2.5,  the  galaxy  emission  should  be  down  by  a  factor  of  5.7.   If 
one  were  to  use  a  receiver  of  comparable  noise  figure  to  the  one  used 
in  this  experiment,  the  galaxy  would  be  well  down  in  the  noise,  while 
the  power  from  the  microwave  background,  having  a  blackbody  spectrum, 
would  be  roughly  constant.   If  there  were  still  no  evidence  of  a 
cosmological  polarization,  the  improvement  in  the  limit  would  be 
roughly  a  factor  of  two. 

Another  major  improvement  could  be  made  by  using  the  second  part 
of  the  dual  mode  transducer  on  the  switch  in  a  two  receiver  system. 
This  would  be  a  form  of  Graham's  receiver  which  would  decrease  the 
effective  system  noise  by  the  square  root  of  two.  Because  the  signal 
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in  the  two  arms  of  the  transducer  are  out  of  phase  at  the  switching 
frequency,  while  the  switch  offset  is  in  phase,  it  also  offers  the 
opportunity  of  balancing  out  the  switch  offset  before  entering  the 
lock-in.   In  this  way  the  immunity  to  gain  drifts  is  enhanced. 

Taken  together,  these  changes  would  give  an  improvement  of  2.5 
in  overall  sensitivity  to  the  background  polarization. 
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APPENDIX  A 

Al.   TIME  DEPENDENCE  OF  THE  SCALE  FACTOES  IN  AN  AXIALLY  SYMMETRIC 
EUCLIDEAN  UNIVERSE  WITH  NO  MAGNETIC  FIELD 
In  order  to  obtain  the  time  dependence  of  the  scale  factors,  one 
has  to  obtain  the  solutions  to  Einstein's  equations  for  the  metric  in 
question,  given  an  appropriate  stress  tensor^  In  this  work  I  deal 

with  an  axially  symmetric  Euclidean  metric  given  by 

2     2       2    2     2        2   2 
ds  =  dt   -  A(t)   (dx  +  dy  )  -  W(t)   dz  .   Metrics  of  this  type 

have  been  considered  by  many  authors,  including  Thorne  (1967) , 

who  has  given  the  following  results  and  whose  notation  I  follow. 

In  the  absence  of  magnetic  fields,  I  assume  a  stress  energy 

tensor  of  the  form, 


Tp     = 

V 

-  P 

0 

0 

0     -1 

0 

/  2 
-p/c 

0 

0 

0 

0 

/  2 
-p/c 

0 

-0 

0 

0 

-p/c  J 

(Al.l) 


for  use  in  the  equations.  The  equations  for  a  quite  general  metric 

57 
have  been  written  down  in  detail  by  Dingle  (1933).    In  particular, 

for  the  above  metric  we  obtain 


8rrGP  _  A    W    AW 

k         A     +   W  +  AW  : 


8rrGP   2A    A2 

T  =  T     +   72 
c  A 


—  ,   (2  equations) 


8nG 


2AW 
AW 


(A1.2) 
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By  subtracting  the  first  equation  from  the  second  we  find  the  following 
expression, 

A   W  .  1?       AW  _  ,    v 

A  "  W+  A2   AW  -  °  '  (A1'3) 

for  which  W  =  A  is  clearly  a  solution.   In  order  to  find  others,  we  will 

make  the  following  two  assumptions.   Since  in  the  limit  of  A  =  W,  we  have 

p/o 

an  Einstein-de  Sitter  universe  for  which  A  =  A~(t/t  )  ■  ,  and  since  the 

case  we  are  considering  is  a  small  simple  departure  from  the  symmetric 

case,  I  will  assume  that  A  is  a  power  law  in  t,  i.e.  A  =  A  (t/t  ) . 

N 

In  addition,  I  will  assume  that  W  is  of  the  form  W  =  A  +  qA  , where 

q  is  some  arbitrary  constant.   In  effect  this  says  that  for  small 
asymmetry, we  expect  a  solution. close  to  that  of  the  Einstein-de  Sitter 
universe  and  for  that  reason  I  ciioose  to  expand  A  in  terms  of  the 
Einstein-de  Sitter  solution.  The  choice  of  q  for  a  constant  is 
unfortunate  and  shouldn't  be  confused  with  q_,  the  deceleration 

parameter.   I  use  it  to  provide  continuity  with  Thome's  results. 

N 

Inserting  A  +  qA  for  W  into  Equation  A1.3  we  find 

q[(l-W)  AN_1  A  +  (1-N2)  AW"2  A2]  =  0.     ,  (Al.U) 

Inserting  the  explicit  form  for  A  gives 

A   \K 

*  [a(*-l)t*N-2  +  (l+N)c2taW-2]  =  0,  (A1.5) 


1  o 
which  reduces  to 

(N+2)o<  =  1.  (A1.6) 
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Next,  taking  the  third  of  Equations  A1.2,  we  note  that  the  left 
side  depends  on  the  density,  which  is  inversely  proportional  to  the 
volume,  meaning  that  Aw  is  a  constant.  Multiplying  this  equation 
through  by  Aw  thus  gives 

A2W  +  2AAW  =  constant,  (A1.7) 

which  becomes,  after  substituting  for  W: 

\ 

'  ?    i\r*p 
3AA  +  qA  A   (1+2N)  =  constant.  (A1.8) 

Again  substituting  for  A  and  using  the  result  A1.6,  we  obtain  the 

expression  for  a, 

A      1 
3t3a-2  _  (3^2)  q  {_o_jj  -   -3  t-2  =  constant?  (A1#9) 

which  is  constant  if  3^-2  =  0  or  a  =   2/3.   This,  with  equation  A1.6, 
implies  that  N  =  l/2,  giving  the  following  forms  for  A  and  W: 

A  =   A  (t/t //3 

°  (ALIO) 

-l/2 
W  =  A  +  qA  '         . 

From  the  above  we  can  now  produce  an  explicit  form  for  the 
differential  Hubble 's  constant: 

Ah  =  -  -  -  .  (Al.ll) 

By  plugging  in  the  appropriate  expressions  for  A  and  W,  we  find 

aA3/2 
Ah  =        °  ,  (A1.12) 

3t  pS] 

o 

2 
which  gives  us  the  result  that  A  WAh  =  Ah  is  a  constant. 


91 


A2.   THE  ANGULAR  DISTRIBUTION  OF  THE  TEMPERATURE  OF  RADIATION  IN  AN 
ANISOTROPIC  EUCLIDEAN  COSMOLOGICAL  MODEL 
This  appendix  will  demonstrate  how,  in  the  absence  of  interactions, 
anisotropically  expanding  Euclidean  models  give  rise  to  an  anisotropy 

in  the  temperature  of  the  microwave  background.   By  Euclidean  models 

2     2       2  2       2  2 
we  mean  models  with  a  metric  of  the  form  ds  =  dt  -  A(t)  dx  -  B(t)  dy 

W(t)2dz2. 

In  order  to  effect  this,  it  is  convenient  to  first  show  that 
p.  -  U.  =  dx./dT  is  conserved  over  the  trajectory  of  each  individual 

blackbody  photon.   In  the  standard  way,  I  define  the  metric  tensor 

iv 

as  g   ,  in  terms  of  which  the  affine  connection  is  given  by 

Ta         =   1/2^  [g   ,  '  +  g  ,   -  g   ,   ]  ,  (A2.1) 

BY  PH  Y    m-Y  P    YB  |i 

in  terms  of  which  the  geodesic  equations  of  motion  of  a  photon  may  be 
expressed  as 

i-  -  ^  **>  •  c*.« 

From  these  definitions,  the  evolution  of  p.  along  the  trajectory 
of  a  photon  may  be  obtained.  Consider  first 


dp. 

dT 

du. 

i 

~~  dT 

d      /          Y\ 
"  dt      (giYU   } 

"Y  h  <eiY) +  elY  if    •  ■  <*2-3) 
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Substituting  from  Equations  A2.2  we  get 
dP 


1 


Y„P  .  „        r_  rY      „«,  B),  (A2.4) 


.  quV  +  g.Y    (-  r   apu  u 


dT     ~  TLY  *  B 
which  by  the  definition  of  the  affine  connection,   Equation  A2.1  becomes 

dP 

dT-=^.Y'pUU      ^Y      (g      }    C^>  p    +   *p'   Y    "   Sjto'  «]  UU 

=  %Y'  puu  "  f  %i»  puu  "  5  %.$>aUUt  *  8^,^  up.  (A2.5) 

In  the  above,  y  is  only  a  dummy  index  and  may  be  replaced  with  a,   giving 

dP. 

i  /  l  saBi  Of  B     i  oB 

dT     =   (%a  '  8    "   2  %d  '  B  )   u  u  "   2  %B''  aa  u  +   2  %p  '  iU  U 

=  T.    u°uP.  (A2.6) 

auB 

Because  dp./dT  has  to  be  symmetric  in  a  and  B  ,  which  are  just 
dummy  indices,  it  must  be  equal  to  the  same  expression  as  follows: 


1    1  rr       „    -1   a  8 

1 —  =  2  [r,  +  r   J  u  uK 

dT    *   caB    Bia 

=  \  [|  &_,,   R  +  g_,    .    -  gifta)   +  I  (g.    ,        +  g      ,    .    -  g.    ,      )3uV 
ia     p         pa     1  18    a        aB    i        is     p 


1    /  \     a  B 

*  (V  i +  V  i}  uu 


=  i'g^,  J  uV   .  (A2.7) 

The  last  line  is  the  needed  result.   It  states  that  if  the  metric 
tensor  is  not  a  function  of  x  ,  then  p.  =  U.  is  a  conserved  quantity 
along  the  trajectory  of  a  photon,  in  the  absence  of  scattering.   In  this 
work  I  deal  with  metric  tensors  which  are  functions  of  coordinate  time 
alone,  which  implies  that  the  photon  momenta  p. ,  i  =  1,  2,  3,  are 
constants  of  the  motion.  This  key  result  allows  us  to  compute  the 
redshift  of  a  photon  emitted  at  an  epoch  t  and  observed  at  the  present 
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time  t  from  the  direction  'n\   Starting  with  the  original  energy  of  the 
photon,  as  measured  by  an  observer  at  rest  with  respect  to  the  x,  y.  and 


z  coordinates,  we  can  write 
(hv)2  =  E2 


=    P      ft*1?      +    P      RJ&   P      +    P    RZi    P 

x  &       SL         y  fe         JL         zg       *l 

~  CT2777  (px)2  +  12777  (V2  +  727"  (pj2] 


A*(t) 


B"(t) 


CT(t) 


(A2.8) 


which  by  conservation  of  the  momenta  p.   is  equal  to 


(hv)2  -  -  [-gi-  (P  )2  +-±~     (P  0)2  +  -|—  (P70)2]. 
A^(t)    X      B^(t)    y      cT(t)    2 


(A2.9) 


For  each  component  of  the  momentum  we  can  show  the  following: 

(pxo)2  -  pxo  er?l  p, 

O 

=  PoPX 
xu 

=  p2  (^  ■  ^)2 

o    v 
=   (hvQ)    (^  •  ^)2, 

which  taken  together  with  Equation  A2.8  gives 


(hv)' 


f4S  .  <o\2-     ° 


ffr  •  *XT+ 


B2(t) 


CnN  •  y)  + 


(A2.10) 


c2(t) 


fc   •  ^)2 


Afc(t) 

By  rearranging  terms  and  taking  square  roots  we  arrive  at  the  following 
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equation  for  the  redshift  z: 


(hv  )' 
v   o 


s* 


1+Z  "  V 


B 


2_  (^  •  <£)2  +    °_  (fr  •  ^)2  +    °        (fr.   ^) 


A2(t) 


1-i 

2 


B2(t) 


c2(t) 


(A2.12) 

With  this  result,  the  temperature  of  the  microwave  background  as  a 
function  of  the  direction  of  observation  can  finally  be  computed.  The 
argument  used  is  similar  to  that  given  for  the  isotropic  case  by 

Tolman  and  rests  on  the  proof  of  Liouville's  theorem  for  general 
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relativity.    Liouville's  theorem  states  that  in  the  absence  of 

interactions,  the  density  of  states,  TJ,  along  the  geodesic  of  a  particle, 

is  a  constant.  Assuming  that  we  are  considering  a  universe  containing 

radiation  with  a  blackbody  distribution  at  recombination,  or  at  least 

at  the  epoch  of  the  last  scattering,  implies 


N  = 


hv/kT  " 
i     '      -1 


=  N  . 


(A2.13) 


using  the  Planck  distribution  characteristic  of  a  blackbody. 
Equality  here  implies  equality  of  the  arguments  of  the  exponentials 
giving 


v 
v    o 

T  "  T 


(A2.3A) 


or 


1+2 


(A2.15) 
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Combining  this  with  Equation  A2.ll,  one  obtains  the  final 

expression  for  the  angular  distribution  of  the  temperature  of  the 

microwave  background  in  an  anisotropically  expanding  Euclidean  universe, 

2  2  2 

m  A  p        p       D  _         p       C  p    ^_1_ 

T  =  - =  T  L—tt-     sin  6  cos  8  +  -=-  sin  9  sin  6  +  -77-  cos  9]~2  . 

0   1+Z      A2  B2  c2 

(A2.16) 

The  result  for  an  axially  symmetric  universe»is  a  specialization  with 

A  =  B  and  C  =  W2  2 

A       ?    W      ?   i 

T  =  T  [-—  sin  9  +  -~-  cos  9]  2  ,  (A2.17) 

°      A  W 

manipulation  of  which  gives 

T  =  T  |  [1  +  -|— ^  -  1  sin29]"2  .  (A2.18) 

o     A  r 

o 

In  the  event  that  the  anisotropy  in  the  expansion  is  small  or,  in  other 


V/A2" 


words,  for  A  W  /AW   =  1.  the  binomial  theorem  allows  us  to  write 
'       o   '    o     ' 

T  =  T   (1  +  e  sin29),  (A2.19) 

o    av  v  ' '  \  si 

where  T   =  TW/W  and  e  =  h   (l  -  A  W  /AT-f  ),  a  small  number  . 
av     '  o        2  v     o   '    o 

From  Appendix  Al  we  can  write  the  explicit  expressions  for  A  and 

W 

A  =  A   (t/t  )2/3 

°     °x  (A2.20) 

W  -  A  +  q  A"2   . 
Insertion  of  these  into  Equation  A2.17  allows  us  to  compare  the  extent 
of  the  expansion  asymmetry  with  the  asymmetry  in  the  temperature 
distribution, 
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e  = 


2A~ 


q  (1  -  (^/t)3)  +  2Aq3/2  (1  -  (Vt)3/2) 


~W 


+  ,2/A  3 


.  (A2.21) 


1  +  2q/Ac 

From  this  it  is  clear  that  the  requirement  that   e  be  small,  so  that 
the  binomial  expansion  in  Equation  A2.l6  is  valid,  is  equivalent  to 
the  parameter  q  being  a  small  number  and  thus,  the  asymmetry  in  the 
expansion  being  small,  which  is  exactly  our  above  assumption. 


A3.   EFFECT  OF  A  SINGLE  THOMSON  SCATTERING  ON  THE  POLARIZATION  OF  THE 
MICROWAVE  BACKGROUND  IN  AN  AXIALLY  SYMMETRIC  UNIVERSE 
When  the  energy  of  a  photon  is  much  less  than  the  rest  energy  of 
an  electron,  its  primary  interaction  with  ionized  matter  is  through 
the  mechanism  of  Thomson  scattering.  This  elastic  form  of  scattering 
can  be  looked  on  as  radiation  by  an  electron  which  has  been  accelerated 
by  an  incident  electromagnetic  wave.   Since  this  is  electric  dipole 
type  radiation,  the  polarization  of  the  outgoing  wave  can  be  derived 
from  the  solution  for  a  simple  radiating  dipole.  If  the  acceleration 
of  the  electron  were  considered  to  be  in  the  a  direction,  the  polar- 
ization of  the  outgoing  wave  would  be  perpendicular  to  the  direction 
of  propagation  and  would  lie  in  the  plane  containing  the  direction  of 
propagation  and  the  direction  of  the  acceleration  of  the  charge  (see 
Figure  A3.l).   To  obtain  E,  the  unit  vector  in  this  direction,  we 
make  the  following  definitions: 

^  -   (sin9  ,  0,  cos9  ) ,  the  line  of  sight  to  the  observer 
a  =  (sin9  cos  $,    sin9sin§,  cos9),  the  direction  of  the 
acceleratd  on  of  the  charge . 

1/   =  the  angle  between  /P  and  /a\ 
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DIRECTION    OF    ACCELERATION 
OF    CHARGE 

LINE    OF    SIGHT    TO    OBSERVER 

POLARIZATION     OF    SCATTERED 
RADIATION 

Z      W    AXIS    OF   THE    UNIVERSE 


FIG.  A3.< 


DIRECTION  OF  SCATTERED  RADIATION  FOR 
THOMSON  SCATTERING 
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From  the  description  given  above  of  the  constraints  imposed  on 
the  polarization  of  the  scattered  radiation,  we  can  write  down  three 
equations  for  the  three  components  of  S.   Perpendicularity  to  the  direction 
of  propagation  r  requires  that  E  •  r  =  0,  while  lying  in  the  plane  con- 
taining r  and  'a" means  that  E  •  (r  x  a)  =  0.   Finally,  being  of  unit 

2    2     2 
length  implies  E  +  E   +  E   =1. 
&     *  x    y     z 

Writing  these  out  in  detail,  we  obtain  \ 


E  •  r  =  E  sinfi,  +  E  cos0  =  0 


x 1    z     1 


E    *    (r  x  a)    =   -cos9_,    sin9   sin$  E     +    (sine,  cos9-cos9    sin9   cos?)   E 
v  1  x        v  1  y 

+   sin9n    sin9   sin§  E    , 
1  z' 

2  2  2 

E+E       +E       =1. 

7  Z  (A3.1) 

Another  expression  which  will  prove  useful  is 

sin  \Jf  =    |r  x  a  |     =  sin  9   sin   §  +    (sin9     cos9   -   cos9     sin9   cos?-)    . 

(A3. 2) 
Solution  of  the  above  equations  yields  the   following  for  the 
components   of  E 

„       cos   9    (sin9     cos9   -  cos9     sin9   cos*) 

E       =  =■ = , 

x 

.2. 
sin  f 

.  2  .2.      .    2./   .    2, 

E       -  sin  9  sm  i/sm  y, 

2 
sin  9   (sin9  cos9  -  cos9   sin9  cos$) 

E  2  =  ± -± g ± ,  (A3. 3) 

sin  i|f 

which  can  now  be  used  to  describe  the  polarization  of  the  scattered  wave, 

In  order  to  put  this  information  in  a  more  useful  form,  I  shall  follow 

Rees '  notation  and  define  the  polarization  of  the  outgoing  wave  in  terms 
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of  the  plane  containing  the  direction  of  propagation,  r  and  the  symmetry 

A 

axis  of  the  universe,  w.  The  direction  perpendicular  to  this  plane 

will  be  called  the  'a'  direction,  while  vectors  lying  in  this  plane 
and  perpendicular  to  r  will  be  said  to  be  in  the  'w'  or  coplanar 
direction.  Applying  these  definitions  to  Equations  A3. 3,  we  find 


2     2 

E   =  E  , 

a    y  ' 

2  2  2 

=   sin   Bsin   §/sin    f    , 

2  2  2 

E=E+E 
W  X  z 

2    2 
=  (sin8  cos8  -  cos8  sin9  cos§)  /sin  ijj  .  (A3.*0 

If  the  incident  radiation  is  polarized  in  the  'a'  direction  and 
has  direction  of  incidence  8  ,  ??,  the  acceleration  of  the  charge  will 
lie  in  the  x,  y  plane  and  be  perpendicular  to  the  direction  of 
incidence.   Thus  8  =  tt/2  and  $  =  8  ±  rr/2.   If  it  is  in  the  'w' 
direction  8  =  8  ±  rr/2  and  $  =  §  +  tt.   Taking  these  cases  separately, 
the  two  components  of  scattered  radiation  (as  seen  by  an  observer 
whose  line  of  sight  r,  makes  an  angle  8,  with  respect  to  the  axis  of 
the  universe  W)  will  be  computed. 

In  Appendix  A2  it  was  shown  that  in  an  anisotropically  expanding 
axial  universe  with  a  small  anisotropy,  the  temperature  of  the 
radiation  can  be  characterized  as  having  a  distribution  of  the  form 
T  =  T   (l  +  e  sin  8),  where  e  is  an  asymmetry  parameter  and  8  is 
the  angle  between  the  direction  of  observation  and  the  axis  W.   In 
this  case,  in  order  to  deal  with  the  polarization,  define  two 
asymmetry  parameters  e  and  e  and  specify  the  temperature  of  the 
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radiation  in  a  given  polarization  in  terms  of  them  by  writing 

2  2 

T  =  T(l  +  e  sin  9_ )    and  T  =  T  (l  +  e  sin  8.,  )  .   Any  change  in  the 
a    v     a     1'  w     v     w     1'  J  ° 

polarization  due  to  scattering  or  other  effects  can  now  be  characterized 

by  its  effect  on  these  parameters. 

The  amount  of  radiation  scattered  in  any  given  direction  is  given 

by  the  product  of  the  incident  flux  and  —  -  =  -^-   amsin  is.   the 

dQ   on  i 

differential  Thomson  scattering  cross  section,  where  a  is  the  total 
scattering  cross  section  and  f  is  the  angle  between  the  direction  of 

acceleration  of  the  charge  and  the  direction  of  scattering.  A  fraction, 

2 
E  ,  of  this  is  scattered  with  polarization  in  the  'a'  direction  and  a 

2 
fraction,  E  ,  with  polarization  in  the  'w ■  direction.  Taking  the  'a' 

direction  first,  the  amount  scattered  into  this  polarization  from 

incident  radiation,  initially  polarized  in  the  same  direction  is 

given  by 

da  „  2 


T   =  f  T  2£  E  <   da. 
a,   J   a  au    a    c 


?T  a  2     2 

T  r  ,,      .  2Q  ,   .  2,  sin  9  sin  i  .    Q  .. 

=  — o (1  +  e  sin  9_)  sin  I ~ d  cos0o  d$0 

on   J  '     a     2  .  d ,  d       d 

sin  v 

3T  a„  ?      ? 

=  -%—   J  (1  +  e&sin  9g)  cos  Sg  2d  cosGp  d$2  (A3. 5) 

T°T    =? 
=  "2-   (1   +ea)' 

In  a  similar  way,  radiation  scattered  into  the  'w'  polarization  is 

given  by 
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da  _  2 


T   =   T  -zk  E   dO, 
w,   J   a  dfi  w    2 


1 


3T  ct  (sin9  cos0  -  cos9  sin9  cos§) 

=  -g^ —  J  (1  +  e&  sin  92)  sin  i|r g d  cos92d§2 

sin  i|t 

T  r>  ?         ?      ? 

=  -q (1  +  e  sin  9„)  cos  9.,  sin  9^  dcos9_d§„ 

on   J  '     a     2       1     2      22 

T  CTT     2     3 
=  ■— -i  cos  B1  (  g  +  ea)-  (A3. 6) 

Starting  with  incident  radiation  polarized  in  the  'w'  direction, 
the  corresponding  results  are 

T  '   =  f  T  &  E  2  dCL 
a?    -J  w  dQ  a    2 

3T  o  2     2 

°   T  p  ,,      .  2.  v   .  2,  sin  9  sin  $   ,    a      .. 
=  -q^ —  J  (1  +  ewsin  e2)  sin  *  d  cos02  d$2 

sin  i|f 

T       r>  9  ?  ? 

=  — g I    (l  +   e   sin  9_.)    cos  9„  sin   §_  d  cos9   2d§0 

on   J  v     w     2'  2      2  2 

T  CTT    1     1 

for  radiation  scattered  into  the  'a'  polarization.   For  the  'w' 
polarization, 

T    '      =    [  T     Jg  E2    d  O 
W„  J      w  dQ     v;  2 


2 

3T  ff  ?  (sin9  cos0   -  cos0   sin9cos§) 


5 (1  +   e     sin  8_J   sin  y 

3n         J    v  w  2 


-  dcos02d$2 
sin  f 

3T  CTT  2  2 

=  —75 (l  +    e      sin   9_)(sin9.,  sin9„-   cos9n  cos0„cos$_)   d  cos9_d§0 

on  Jv  w  2/x  1  2  1  2  2 '  22 


^r1    li'^J^ si"\< I  ♦  ii   V> 


(«.8) 
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Adding  these  results  gives  the  total  contribution  to  each 
polarization 

T        =  T        +  T        =  Tam    (1  +   h  e  +     ±-      e    ) 
a  a  a  T    v  2     a       10        w 

T        =  T       +T'=Tarp[(l  +  ^e+Tx      e   )   +    (  :L      e-|e    )sin26,  ]. 
w  w  w  T      v  d     a       10        w'         v  30        w      ■=     ay  1 

(A3. 9) 
If  we  now  re-express   equations  A3. 9  in  terms,  of  new  asymmetry  parameters 

e    '    and  e    '   we  obtain  to  first   order  in  e     and  e    .the  following  solutions 
aw  aw' 

for  the  new  parameters : 

e   *  =  0 

a  (A3. 10) 

e    '    =    I-      e   -     i      e    , 
w  30        w       2        a' 

where  the  new  average  temperature  is  T '  =  To"m(l  +  h   e  +  —      e  ). 

°     *  Tv    2  a    io   w 

In  matrix  form  this  can  be  written 


0     0 

(A3. 11) 


This  last  expression  is  the  final  result  for  the  effect  of  a  single 
Thomson  scattering. 


Ai+.   TERRESTRIAL  OBSERVATION  OF  A  COSMOLOGICAL  LINEAR  POLARIZATION 

In  previous  sections  I  have  shown  how  cosmological  causes  can  give 

rise  to  a  linear  polarization  of  the  microwave  background  of  the  form 

(T  -  T  )  =  (T  -  T  )    sin2e,  (Ak.l) 

v  w   a.'        v  w   a  max      '  v     ' 

where  9  is  the  angle  between  the  line  of  sight  of  the  observer  and 

the  direction  defined  by  the  symmetry  axis  of  the  universe.  Now,  the 

question  is  what  form  does  this  polarization^  take  for  a  ground  observer 

trying  to  measure  its  effect?  In  my  case  I  will  assume  the  following 

measurement  scheme.   A  microwave  polarimeter  will  be  aimed  at  the 

zenith  and,  as  the  earth  rotates,  will  be  allowed  to  sweep  out  a  circle 

of  constant  declination  in  the  sky.  Given  this  situation,  we  need  to 

find  the  variations  in  the  signal  produced  by  the  rotation  of  the 

earth.   We  can  diagram  the  experimental  setup  as  in  Figure  A^.l, 

where  E  and  E  ,  the  unit  vectors  in  the  perpendicular  and  coplanar 

directions  as  defined  in  Appendix  A3,  have  been  drawn  in. 

By  stating  their  relationships  to  the  directions  R  and  W  and  to 

the  normal  to  the  R,  W  plane  R  x  W,  we  can  develop  equations  for  the 


determination  of  E  and  E  .   First,  because  both  of  .them  are  parallel 

aw  * 

to  the  polarization  of  incident  radiation,  they  must  be  perpendicular 
to  the  line  of  sight,  R.  Next,  E  lies  in  the  S,  W  plane  and  must  be 
perpendicular  to  its  normal  while  E  is  perpendicular  to  this  plane 
and  thus  must  be  perpendicular  to  W.   Finally  each  of  them  must  be  a 
unit  vector.   In  vector  notation  these  conditions  become 
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N- NORTH  CELESTIAL  POLE 

W-  AXIS  OF  THE  UNIVERSE 

R-  LINE  OF  SIGHT  OF   THE  OBSERVER 

8  =  7T/2-S 

a  =  cl0-t  (where  r  is  local  sidereal  time 

in  angular  units) 


FIG.  A  4.1 

RELATION  BETWEEN  THE  AXIS  OF  THE 
UNIVERSE,  THE  ROTATION  AXIS  OF  THE 
EARTH  AND  THE  LINE  OF  SIGHT  TO  THE 
OBSERVER  • 
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E      '    §  =  0 
w 


E      •    g  =  0 
a 


E      •    (W  x  R)    =  0 


w 


2  2  2 

E+E+E  =1 

wx  wy  wz 


E      •   W  =  0 
a 

2  2  2 

E+E  +E  =1, 

ax  ay  az 


(AU.2) 


while  written  out  in  detail,  they  are 


E    sin  6  +  E   cos  6=0 
wx  wz 


E   sin6  +  E  cos6  =  0 
ax        az 


E    sin  §  sina  cos6  +  E  (cos9sin6 
wx  wy 

-sinScosacosS)  +  E  sin§  sino-  sin6  =  0 
wz 

2      2      2 
E    +  E    +  E    =1 
wx     wy     wz 


Taking  these  with  the  expression 


E     sin$  cosa  +  E     sin9    sine* 
ax  ay 

+  E     cos$  =  0 
az 

2  2  2 

E+E  +   E  =  1. 

ax  ay  az 


(AU.3) 


p  p p 

|w  x  R|   =  sin  9   =  sin$  sin  a  +    (cos§  sin6    -  sin$  cosa  cos6)    ,  (Ak.k) 

leads  to  the  following  solutions   for  the  unit  vectors: 


E 


wx 


cos6  (cos?  sin6  -  sin$  cosa  cos6) 
sin9 


ax 


cos6    sin§  sincr 
sin9 


„  sm§  sina 

E        =   +   r~5 

wy  sin9 


ay 


sin&  cos  I   -   sin5  coscy  cos 6 
sin9 


sin6    (cos§   sin6    -   sin$  cosa  cos6) 
sin9 


wz 


E 


az 


sin5    sin9    sincr 
sin9 


(A^.5) 

this  allows  us  to  express  the  resultant  field  seen  by  the  observer  as 


E  =  E  E  +  E  E   , 
a  a    w  w  ' 


(AU.6) 


with  the  amplitudes  E  and  E  complex  numbers  and  functions  of  t.   In 
order  to  resolve  these  into  components  relative  to  directions  with 


io6 


respect  to  the  north  celestial  pole,  we  notice  from  Figure  A^+.l  that 
our  choice  of  coordinates  implies  that  the  East-West  component  of 
E  is  the  E  component  while  the  North-South  component  is  given  by  the 
projections  in  the  plane  of  the  other  two  or 


E  =  E  sin  6  -  E  cos  6 
ns   z         x 


E   =  E  . 
ew    y 


(A4.7) 


\ 


By  applying  a  rotation  to  these  two  components,  the  field  can  be  expressed 
relative  to  any  direction  we  choose.   For  an  axis  at  an  angle  (3  with 
respect  to  North,  we  can  give  the  components  with  respect  to  the  new 
axis  as 


cosp  sinp 
-sin(3  cosp 


ns 


(All. 8) 


E 


ew, 


For  a  polarimeter  of  the  type  I  am  using,  the  important  quantity 

is  T   -  T   ,  which  is  equal  to  C ( | E  I  -  |e  |  ).   The  quantity  in  brackets 

can  be  written  in  terms  of  E   and  E   as 

ns      ew 


IE.  I2-    IeJ2  =    (I    E      I2-    |E      I2)    cos2p  +   2    |E      •    E     *|    sin26, 
1   1 '         i    2 '  x '      ns  '         '   ew '  K  '   ns  •      "t-t    ' 


ew 


(Ak.9) 


or  in  terms  of  E     and  E     as 
aw 

l_    i2      .„    ,2    rlEw'    "    lEa'    I,         ,      .    T  .    .         -         r,2        .    2-     .    2— , 

|E    |    -    |E    J   =j L(cos$  sine  -  sini  cosa  cos6)    -   sin   $  sin  a  J 


sin  0 

i            *i 
2  E   -E 
1  w    a  '  /       -     .  t        .    -        —        —«.,.— 
+  p (cos?  sino   -   sin$  cosa  coso)sin$  sma 

sin  6 


cos2p 


+2 


(k, I2-  K,l2) 


(cos$  sin6    -   sin?-  cosa  cos6)    sin*  sine* 


sin 


|E   •    E   * 

1   w        a    '  2  2  -  - 

P         '  l_sin   i  sin  a  -    (cos$  sin6-sin$  cosa  cos6) 
sin  9 


2 

(AU.10) 


sin28. 
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In  a  practical  situation  where  the  signal  is  integrated  for  times  much 
greater  than  the  reciprocal  of  the  receiver  bandwidth,  we  are  taking  a 
time  average  of  the  above  expression.   Because  the  signals  in  the  two 
polarizations,  E  and  E  ,  are  independent,  the  time  average  of  their 
product  is  zero  and  we  obtain 

i         _  _    _  2     2     2— 

T  -  T^  =  (T  -  T  )    <|7cose  sin6  -  sin?  cosa  cosf)  -  sin  $  sin  oil   cos2B 
1   2   v  w   a  max  1^ 

%  (A4.11) 

+  2(cos§  sin6  -  sin*  cosa  cos6)  sin§  sina  sin2(3>  , 

where  we  have  multiplied  through  by  the  constant  C  and  made  the  sub- 

2 

stitution  (T  -  T  )/sin  8  =  (T  -  T  )    .  This  result  is  the  signal 
v  a   wx/        v  a   wmax 

desired.   Replacing  a     and  6  with  their  equivalents  in  terms  of  the 
local  sidereal  time  t   (this  is  equal  to  the  right  ascension  of  the 
zenith),  the  right  ascension  of  the  axis  of  the  universe,  a   ,  and 
the  declination,   6  ,  and  doing  some  manipulation  we  arrive  at  the 
following:  which  is  the  polarization  signal  received  by  a  switched 
polarimeter  pointed  at  the  zenith  which  has  its  positive  axis  making 
an  angle  p  with  respect  to  North.  Here,  positive  is  defined  as 
the  T  direction: 

T  -  T0  =  (T  -  T  )    )[cos26  (1-|  sin2e)  +  sin2$  (l  -  J  cos26)  cos2  (a  -   t) 
1   2   *  w   a  max  J       v    2      '         v    2  x  o 

1  2 

-  77  sin2$  sin26  cos(a  -  t)j  cos2B  -  [sin  §  sinS  sin2(a  -  t) 

2  o  o 


-  sin2§  cos6  sin  (a  -   t)]  sin2B 


o 


(AU.12) 
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APPENDIX  B 

Bl.   THE  ANTENNA 

As  I  had  nothing  to  do  with  design  of  the  antenna,  I  will  not  go 
into  detail  on  the  techniques  of  design  but  will  only  present  a 
description  of  its  operating  characteristics.  Interested  readers  can 
refer  to  the  references,  in  particular  the  paper  by  Roll  and  VJilkinson, 
for  further  information. 

The  antenna  used  was  a  conical,  optimum  gain  horn  with  an  aperture 
choke  for  suppression  of  the  backlobes.  Provision  was  made  as  shown 
in  Figure  Bl.l  for  the  insertion  of  a  calibration  probe  in  the  side  of 
the  horn,  and  this  feature  was  in  place  during  about  two  months  of 
observations.  The  full  width  between  half  power  points  is  about  15 
degrees  which  gives  the  polarimeter  a  resolution  of  one  hour  in  right 
ascension.   Figure  B1.2  is  a  plot  of  antenna  gain  as  a  function  of  angle 
from  the  horn  axis  for  illumination  by  orthogonal  polarizations.  As 
the  procedure  varies  slightly  from  that  used  for  a  power  radiometer, 
a  little  explanation  is  in  order.   In  making  the  measurement,  the 
polarimeter  was  mounted  horizontally  on  a  rotating  table  with  the 
apex  of  the  horn  at  the  center  of  rotation.   One  polarization  sensitive 
axis  of  the  polarimeter  was  aligned  horizontally  and  the  other 
vertically.   Incident  radiation  was  provied  by  a  klystron  radiating 
through  rectangular  waveguide  and  a  standard  gain  horn,  which  gave  a 
linearly  polarized  source.   Horn  to  horn  distance  was  a  little  over 
twenty  meters,  giving  a  good  approximation  to  a  plane  wave  at  3-2 
centimeters.  Two  scans  were  made,  one  with  the  incident  polarization 
horizontal,  which  corresponded  to  the  E-plane  measurement  of  a  single 
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polarization  microwave  horn,  and  the  other  with  polarization  vertical, 
corresponding  to  the  K  plane  measurement.  Because  of  these  analogies  I 
have  used  the  terms  E  and  H  plane  in  labeling  the  resultant  plots.   The 
appearance  of  the  E  plane  measurement  is  very  similar  to  what  one  would 
find  for  a  power  radiometer,  while  the  H  plane  measurement  has  a  striking 
difference.   For  some  angles  the  sign  of  the  signal  changes,  indicating 
a  coupling  between  the  polarizations.   It  is  not  known  whether  the 
observed  effect  is  due  to  external  reflections  or  to  coupling  within 
the  polarimeter  itself.  Because  the  effect  occurs  at  levels  of  -35db 
and  below,  it  has  little  effect  on  the  experiment  and  was  thus  not 
pursued  further . 

The  primary  purpose  of  the  measurements  was  to  produce  Figure  B1.3, 
which  is  the  difference  between  the  E  and  H  gain  patterns.  This 
difference  gives  directly  the  response  to  unpolarized  radiation,  which 
can  then  be  used  to  compute  the  expected  contribution  from  ground 
radiation.  The  antenna  temperature  due  to  any  source  falling  within 
the  beam  pattern  is 

T  =  J  TP(fi)dcyJp(fi)dO,  •  (Bl.l) 

where  T  is  the  equivalent  blackbody  temperature  of  the  source  at  3.2 
centimeters  and  P(Q)  is  the  gain  of  the  antenna.  The  denominator  can 
be  evaluated  by  noting  that  calibration  is  performed  by  covering  the 
entire  antenna  pattern  with  a  source  at  temperature  T  and  calling  the 
resultant  antenna  temperature  T.   For  a  rough  approximation  of  ground 
radiation  we  can  use  the  form 
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Ta=  TPAV(n)An/Jp(n)dn.  (el. 2) 

At  the  observation  site  on  the  roof  of  Jadwin  Physics  building, 
the  most  prominent  terrain  feature  and  the  expected  source  of  most  of 
the  ground  radiation  is  a  structure  on  the  roof  adjacent  to  the  polarimeter. 
This  structure  houses  the  elevator  motors  and  is  approximately  two 
stories  high  or  an  angle  of  about  a  half  radian  above  the  horizontal. 
Its  width  is  about  a  third  radian  ,   giving  it  a  size  of  about  one 
sixth   steradian  located  between  an  angle  of  sixty  and  ninety  degrees 
from  the  horn  axis.   In  this  position  the  antenna  gain  from  Figure  B1.3 
is  about  -37db  for  unpolarized  radiation.   If  we  take  29O  K  as  an 
average  temperature  for  the  structure  we  obtain 

TP   (n)An        a  „  ,9Qr0 

as  an  appropriate  contribution  to  the  signal  from  ground  radiation  for 
an  unshielded  antenna. 

During  the  experiment  the  antenna  was  shielded  from  ground  radiation 
by  a  shield  which  extended  from  below  a  zenith  angle  of  ninety  degrees 
to  a  zenith  angle  of  thirty  degrees.  Measurements  of  its  effectiveness 
using  the  klystron  source  were  done  for  zenith  angles  of  between  eighty 
and  ninety  degrees.  The  average  reduction  in  sensitivity  produced  by 
the  shield  in  this  range  was  12.3db.   If  we  take  this  figure  and  apply 
it  to  the  above  calculation,  we  find,  assuming  that  the  shielding  factor 
is  valid  over  this  range,  that  a  reasonable  figure  for  the  expected 
contamination  from  ground  radiation  is  7-9  radeg.  K.   In  the  first  two 
months  of  the  experiment  I  noticed  a  DC  level  of  polarization  of  6.66 
mdeg.  K.,  which  agrees  with  the  rough  estimate  above.   I  felt  that  this 
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was  probably  not  unreasonable.   An  order  of  magnitude  measurement,  using 
the  shielding  factor  and  a  piece  of  eccosorb  as  a  source  inside  the 
shield,  showed  rough  agreement.   In  fact,  this  measurement  was  correct 
at  that  time.   It  wasn't  until  I  later  decided  that  I  would  not  need 
the  calibration  probe  and  therefore  removed  it  from  the  horn  that  I 
discovered  that  the  high  degree  of  sensitivity  to  ground  radiation  was 
due  to  the  perturbing  influence  of  the  probe. <.  In  the  months  that 
followed  the  level  of  DC  polarization  dropped  to  less  than  1.5  mdeg  K. 
By  scrutinizing  the  level  as  a  function  of  time  I  was  able  to  localize 
the  effect  to  the  approximate  time  of  the  probe  removal.   Subtraction 
of  this  effect  is  covered  under  data  analysis. 

From  the  above  discussion  one  can  see  that  effects  from  ground 
radiation  are  no  trivial  matter  and,  in  fact,  form  part  of  the  limit 
on  the  sensitivity  of  this  whole  experiment.  This  was  one  of  the  key 
reasons  for  observing  only  the  zenith  where  these  effects  would  be 
minimized. 
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B2.   THE  SWITCH 

In  order  to  study  the  effect  of  the  switch  in  detail,  I'll  first 
consider  it  as  an  ideal  Faraday  rotator.  This  is  done  in  order  to 
determine  precisely  the  signal  as  a  function  of  the  direction  of  the 
incident  polarized  radiation  and  of  the  magnetic  field  used  to  drive 
the  switch.  A  convenient  reference  system  for  specifying  the  direction 
of  polarization  of  incoming  radiation  is  the  symmetry  plane,  which 
contains  the  axis  of  the  ferrite  and  bisects  the  right  angle  formed 
by  the  orthogonal  directions  defined  by  the  dual  mole  transducer  (See 
Figure  B2.l).   Incoming  radiation  will  be  specified  in  terms  of  the 
magnitude  of  its  electric  vector  and  angle  between  this  vector  and  the 
symmetry  plane  9  . 

After  propagating  along  the  ferrite  the  E  vector  has  been  rotated 
by  an  amount  determined  by  the  longitudinal  magnetic  field  applied  to 
the  ferrite.  In  our  case  we  are  driving  the  switch  with  a  sinusoidally 
varying  current  and  thus  the  angle  will  be  equal  to  Acosat  +  B  where 
A  is  the  maximum  rotation  angle  and  B  is  the  rotation  caused  by  any 
DC  magnetic  fields  present.  Combining  these  we  get  an  expression  for 

— * 

the  angle  between  the  direction  of  the  E  vector  and  symmetry  plane  of 
the  switch  just  before  the  radiation  enters  the  dual  mode  transducer 

9  =  9  +  A  cosuut  +  B.  (B2.1) 

o  v     ' 

Looking  at  the  signal  from  one  arm  of  the  dual  mode  transducer,  we 

obtain  the  following  for  the  electric  field  in  the  waveguide  entering 

the  receiver 

Earm2=EcoS  (J"9)'  <B2'2) 
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In  the  detection  process  the  receiver  produces  a  signal  proportional 
to  the  square  of  the  electric  field,  which,  in  turn,  is  proportional 
to  the  temperature  of  the  received  radiation.   In  this  way  the  signal 
leaving  the  receiver  can  be  expressed  in  terms  of  the  temperature  of 
the  incident  radiation  as 

2  ,TT 


S  a   T  cos 


(it  "  6).  (B2.3) 


Substituting  for  9  and  expanding,  we  obtain  * 

CT    ( 

S  =  ~   jl  +  cos2B  [sin29   cos(2A  cos  (cut))  +  cos26   sin(2A  cos(uvt))] 

+  sin2B  [cos29   cos(2A  cos  (tut))  -  sin29   sin(2A  cos(tut))]( 

(BZ.k) 
for  the  signal  present  at  the  input  to  the  lock-in  amplifier, 
where  C  is  a  proportionality  constant  depending  on  the  gain  of  the 
receiver. 

The  action  of  the  lock-in  is  to  extract  the  Fourier  component  of 
the  fundamental  at  the  switching  frequency.  The  labor  can  be  eased 
considerably  by  noting  that  most  of  the  terms  in  the  above  expression 
do  not  have  this  period.   The  first  term  is  DC  and,  of  course,  gives 
no  contribution.  The  terms  in  cos  (2Acos(tut))  are  symmetric  about 
tut  =  0  and  have  a  frequency  twice  that  of  the  driving  current.  For 
A  =  n/k   the  waveform  is  similar  to  Figure  2.6b,  which  also  does  not  have 
a  Fourier  component  at  the  switching  frequency.  All  that  remains  are 
terms  in  sin(2Acos(tut)),  which  give  signals  like  Figure  2.6a  and  are 
the  only  terms  that  contribute.   By  adjusting  the  phase  of  the  lock-in 
for  maximum  signal,  we  account  for  phase  shifts  in  the  apparatus  and 
find  for  the  Fourier  component  and  thus  the  output  of  the  lock-in 
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V,      ,     .   =  ^rr-  [cos2B  cos29    -   sin2B  sin20    ]    (-)     \"     sin(2A  cosuut)    coscutdt. 

lock-in   2  o  o  ^n'  Joo  '  ' 

(B2.5) 

where  D  is  a  constant  depending  on  the  gain  of  the  lock-in. 

The  solution  to  this  integral  is  the  Bessel  function  J  (2A)  which  gives 

CUT 
V,   .  .   =  ~  cos2f9  +  B)  2J_ (2A).  (B2.6) 

lock-in    2      '  o       lv  v    ' 

In  order  to  maximize  the  output,  J  (2A)  has  to  be  maximized.   Setting 
dj  (2A)/dA  =  0,  we  note  that  the  first  zero  of  J  ' (x)  occurs  at  x  =  1.841 
and  corresponds  to  a  maximum.  This  implies  that  for  the  greatest  signal, 
the  maximum  rotation  A  is  52.7  degrees.  The  value  for  J  (2A)  is  . 58187, 
or  in  terms  of  the  output  voltage  , 

Vlock-in  =  lT  C1"16^  -s2(eo+B),  (B2.7) 

which  differs  from  that  obtained  for  a  purely  sinusoidal  signal  by  the 
factor  1.164.  This  has  to  be  included  when  interpreting  the  results. 
The  minimum  detectable  signal  for  the  polarimeter  is  also  altered  to 

account  for  the  change  in  effective  amplitude  of  the  modulating  signal 

T 

AT  .   =2/2   — ^ .  (B2.8) 

110X11  (l.l64)/AvAt 

Returning  now  to  equation  B2.6,  we  note  that  if  the  residual  fields 

represented  by  the  rotation  B  are  DC  in  nature,  the  only  effect  is  a 

change  in  the  effective  position  of  the  symmetry  plane  of  the  switch, 

an  effect  which  can  easily  be  calibrated  out.   If  there  is  an  AC 

component  due  to  rotation  or  other  movement  of  the  device,  we  can, 

given  that  it  is  small,  give  the  following  expression  for  the  lock-in 

output  to  first  order  in  B  , 
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V,   .  .   =~^  (1.16U)  [cos20'  -  2B'  sin29'],  (B2.9) 

lock-in   2  v  o  o 

■where  0'  =  9  +  B,   and  B'  is  the  AC  component  of  B. 
o   o   dc 

In  my  case,  where  the  radiometer  is  mounted  vertically  and  rotated 
about  the  vertical  axis,  the  maximum  possible  variation  in  the 
longitudinal  magnetic  field  is  that  due  to  changes  in  the  projection 
of  the  horizontal  component  of  the  earth's  field  on  the  axis  of  the 
ferrite,  caused  by  deviation  in  the  alignments  of  this  axis  with  the 
vertical.  The  horizontal  component  of  the  earth's  field  in  Princeton 
is  about  .2  Gauss.  Because  a  misalignment  would  have  a  small  projection 
in  the  horizontal  direction,  .2  Gauss  is  clearly  an  upper  limit  on  this 
effect  for  any  orientation.   Therefore,  an  upper  limit  on  the  change 
is  twice  this  number.  My  best  estimate  of  the  driving  field  of  the 
switch  was  5.3  Gauss  and  came  from  wrapping  a  coil  around  the  switch 
throat  and  measuring  the  induced  voltage.  This  number  should  be  less 
than  the  maximum  driving  field  because  the  coil  was  about  five  centimeters 
above  the  center  of  the  ferrite,  thus  it  should  provide  a  conservative 
estimate.   If  one  makes  the  assumption  that  the  Faraday  rotation  in  the 
ferrite  is  a  linear  function  of  field  strength,  then,  the  earth's  field 
can  amount  to  no  more  than  a  7.5fa  effect.   In  actual  practice,  magnetic 
shielding  of  the  svtitch  reduced  this  effect  by  another  factor  of  150 
and,  hence,  to  insignificance. 

Another  and  more  important  effect  and  the  one  which  actually  made 
the  shielding  necessary  was  the  change  in  the  absorption  coefficient  of  the 
ferrite  with  longitudinal  magnetization.  Figure  B2.2a  is  an  abosrption 
curve  for  a  typical  ferrite  used  in  microwave  work.    In  order  to 
demonstrate  the  effect  we  are  considering,  it  is  only  necessary  to 
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note  that,  in  general,  the  curve  is  non-linear  and  in  any  particular 
ferrite  one  is  somewhere  on  the  curve  in  a  non-linear  region.   In  a 
device  like  the  switch,  where  an  alternating  field  is  applied,  the 
absorption  varies  sinusoidally  with  the  applied  field  around  the 
zero  field  value  (Figure  B2.2b).   If  we  change  the  zero  field 
absorption  by  biasing  the  switch  with  an  external  magnetic  field, 
such  as  the  earth's  field,  the  magnitude  of  the  sinusoidal  component 
of  the  absorption  changes  due  to  the  change  in  slope  of  the  curve. 

Because  I  am  detecting  the  signal  with  a  lock-in  amplifier, 
I  am  only  sensitive  to  changes  in  the  absorption  at  the  switching 
frequency.  In  order  to  display  this  explicitly,  I  will  write  the 
absorption,  a,   in  the  following  form  ; 

or  =  a^   (M)  +  aAC    (M)  coswt,  (B2.10) 

where  the  departure  from  exact  sinusoidal  dependence  to  non-linearity 
of  the  absorption  curve  has  been  neglected,  as  it  is  only  the  first 
harmonic  I  am  sensitive  to. 

In  passing  through  the  switch,  unpolarized  radiation  undergoes 
absorption,  which  is  proportional  to  the  producL  of  the  absorption 
coefficient  and  the  temperature  of  the  radiation: 

P-  =  -  a  T,  .  (B2.il) 

dx       in  v     ' 

In  addition,  there  is  a  gain  in  power  due  to  emission  in  the  ferrite, 

giving  rise  to  a  signal  proportional  to  the  emission  coefficient  (3 

times  the  temperature  of  the  ferrite  T  : 

f-  =   BTo.  (B2.12) 


122 


TYPICAL    FERRITE    ABSORPTION  CURVE 


a 


MAGNETIZATION 


FIG.  B2.2  (a) 


a 


2a 


AC  (2) 


AC   ABSORPTION  FOR  TWO  VALUES 
OF  DC  BIAS    FIELDS  82 


2a 


AC(1) 


DRIVING  MAGNETIZATION 
OF  SWITCH 


■*-    MAGNETIZATION 


FIG.  B2.2  (b) 


TYPICAL  FERRITE  ABSORPTION  CURVES 


123 


The  relation  between  <x     and  (3  can  be  established  from  consideration 
of  detailed  balance.  If  the  input  radiation  were  made  to  be  equal  in 
temperature  to  the  switch,  the  expressions  in  B2.ll  and  12  would  have 
to  be  equal  in  order  that  thermal  equilibrium  be  maintained,  which 
would  imply  that  a   =  B.  We  can  write  for  the  change  in  temperature 
of  the  radiation  in  passing  through  the  switch, 

rlT  ^ 

=~-   =  (T  -  T.  )  a,  (B2.13) 

dx   v  o   in 

which  has  as  its  solution, 

T   .=  T.   e"T  +  T   (1  -  e"T),  (B2.lU) 

out   in       o  v        '  v     ' 

where  t  is  a  times  the  length  of  the  ferrite  and  is,  in  effect,  the 

'optical  depth'  of  the  switch.  For  t  small  this  result  can  be  given 

to  first  order  in  t  by 

T   .  =  (T  -  T.  )  t.  (B2.15) 

out   x  o   in'  v 

Under  normal  operating  conditions  T  -  T.   is  about  300  K.  With  the  AC 
*  o   in 

portion  of  the  absorption  this  gives  an  offset  signal  at  the  switching 
frequency  of 

T    =  300  t   (M)  coso)t°K.  (B2.16) 

OUT*  ElC 

Actual  measurement  of  this  quantity  gave  a  value  of  20  to  100  mdeg.  K 
depending  on  the  temperature  of  the  ferrite.  This  evidence  gives  two 
pieces  of  important  information.  First,  it  says  that  t  (m)  for  the 

9.C 

-h 

switch  is  of  the  order  of  3.0  x  10    and  thus,  direct  modulation  of  the 
input  polarized  signal  by  the  absorption  in  the  ferrite  can  be  neglected. 
Secondly  and  more  important,  it  demands  that  one  look  carefully  at  the 
variation  of  magnetization  with  radiometer  orientation  due  to  the 
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earth's  magnetic  field. 

Direct  measurement  of  the  sensitivity  of  the  switch  to  perturbing 
longitudinal  magnetic  fields  showed  a  change  in  switch  effect  of 
156  ±  3  mdeg.  K/Gauss  for  a  change  of  2.1+  Gauss;  and  a  change  in 
offset  of  130  ±5.5  mdeg.  K/Gauss  for  a  change  of  1.1  Gauss. 
Application  of  up  to  11  Gauss  differences  perpendicular  to  the  ferrite 
axis  failed  to  show  any  effect  to  the  same  order  of  accuracy  as  the 
longitudinal  results.  In  order  to  gauge  the  effect  of  these  offsets, 
assume  that  the  ferrite  axis  was  ten  degrees  from  the  vertical.  The 
projection  of  the  horizontal  component  of  the  earth's  field  on  the 
ferrite  would  result  in  a  maximum  change  in  the  longitudinal  field 
of  ,k   sin  (10  )  or  about  .07  Gauss.   From  the  above  data  this  would 
result  in  a  variation  in  the  signal  of  up  to  9-3  mdeg.  K,  a  whopping 
signal!   In  the  early  days  of  debugging  the  apparatus,  signals  of  this 
order  were  actually  observed  and  on  the  basis  of  this,  it  was  decided 
to  shield  the  switch  with  high  permeability  material.  Tests  run  on 
the  actual  shield  used  showed  that  the  shield  reduced  the  variation 
in  the  horizontal  component  of  field  by  a  factor  of  150  and  all  fields 
by  better  than  a  factor  of  100.  In  this  way  the  problem  of  stray  magnetic 
fields  was  also  eliminated. 

The  final  item  to  be  considered  is  the  change  in  switch  offset  with 
temperature  of  the  ferrite.   Good  numbers  on  this  were  available  from 
the  actual  data  taken  in  the  experiment  itself.   In  addition  to  measuring 
polarization  of  the  microwave  background,  I  measured  the  switch 
temperature  as  a  function  of  time.   By  averaging  together  all  the  data 
at  a  given  temperature,  I  can  obtain  the  DC  offset  as  a  function  of 
temperature  (see  Figure  B2.3).   It  may  seem  odd  that  I  was  able  to  observe 
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a  significant  change  in  temperature  in  spite  of  controlling  the 
environment  inside  the  experiment  shed.  As  it  turned  out,  heating  of 
the  horn  antenna  by  the  sun  during  the  day  and  conduction  of  this  heat 
to  the  switch  by  the  metal  waveguide  caused  an  uncontrollable  variation 
in  the  switch  temperature.  Compensating  for  extreme  changes  in 
temperature  during  a  single  day  was  also  hard.   My  controlling  apparatus 
had  a  limited  capacity  and  had  to  be  biased  into  a  proper  range  by 
controlling  the  ventilation  to  compensate  for  the  season  of  the  year. 
Abnormal  fluctuation  caused  the  temperature  to  depart  from  the  controll- 
ing range. 

Luckily,  the  increase  of  offset  with  temperature  variation  is  not 
serious  as  long  as  the  variations  occur  on  the  time  scale  of  a  few 
integration  times.  The  analysis  process  subtracts  out  offsets  and 
under  these  conditions  the  temperature  induced  ones  would  be  eliminated. 
The  worst  conditions  are  heating  by  the  sun  on  a  hot  day  where  the 
temperature  can  go  up  by  as  much  as  3  degrees  C/hr  during  the  heating 
phase  and  drop  by  about  2  degrees  c/hr  in  the  cooling  phase.   I 
ignore  the  heating  because  data  taken  during  this  time  is  contaminated 
by  the  sun  and  not  used  anyway.  The  variation  during  the  cooling  phase 
and  Figure  B2.3  can  be  used  to  estimate  an  offset  change  of  approx- 
imately 15  rndeg  K/hr  during  the  first  couple  of  hours  after  the 
sun  goes  down  behind  the  ground  shield.  This  effect  did  not 
contribute  enough  solar  signal  in  the  analysis  to  warrant  discarding 
additional  data. 
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SPECIFICATIONS  OF  THE  SWITCH 


Ferrotec  Rotational  Circulator  Switch  (with  1  dual  mode  transducer  removed) 


Model 

R-165LS 

Serial  Number 

11 

Insertion  Loss 

.1  db 

Frequency- 

8.2-10.0  GHz 

Current  for  opi 

;imum  switch 

rotation 

6k   ma 

Optimum  switch 

rotation 

52.7° 

Sensitivity  to 

Long. 

130  ±  5.5  md 

Sensitivity  to 

B± 

0  ±  10.0  mde 

Offset  vs.  temp.  +6.7  mdeg/°C  temp.  >  3o°C 

•4.7  mdeg/°C  temp.  <  36°C 
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B3.   THE  RECEIVER 

The  receiver  consisted  of  the  following  components: 
MIXER 

BHG  Electronics  Laboratory  Model  MP8/4AAEFI 

Serial  Number  5-163-1 

Frequency  8.5-9-6  GHz 

IF  Frequency  60  MHz, 

IF  Bandwidth  70  MHz 

Power  Gain  29  db 

Noise  Figure  7.2-7.3  db 

IF  AMPLIFIER 

RHG  Electronics  Laboratory  Model  EBT101MGC 

Serial  Number  5-163-2 

Center  Frequency  60  MHz 

3  db  Bandwidth  48  MHz 

Maximum  Power  Gain  75  db 

LOCAL  OSCILLATOR 

Microwave  Associates  Gunn  Oscillator  Model  MA86101 

Serial  Number  1190 

Frequency  9.37  GHz 

Power  Output  10  mw 

A  useful  comparison  can  be  made  between  the  quoted  noise 

figure  of  the  mixer  given  above  and  the  actual  system  noise. 

A  total  temperature  of  lkk5   K  was  observed  by  covering  the 

polarimeter  with  a  29*+  K  absorber  and  measuring  the  variance 

of  the  data.  This  gives  a  system  temperature  of  l'-iU5-2^U  =  1151  K, 

It  must  be  noted  that  this  number  is  referenced  to  the  290  K 
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temperature  difference  between  the  sky  and  the  piece  of  eccosorb  used 
for  calibration.   Because  of  insertion  loss  in  the  front  end  of  the 
polarimeter,  the  size  of  the  calibration  load  is  reduced  in  magnitude 
by  the  time  it  reaches  the  mixer.  A   rough  insertion  loss  can  be 
obtained  from  summing  the  following  contributions : 


The  Switch 

.1  db 

The  Circulator 

.2  db 

The  Isolator 

.8  db 

Waveguide  joints 

.2  -  .3  db 

Total 

l.k   db 

Reducing  the  input  reference  load  by  this  amount,  we  find  that  the 
noise  figure  is  based  on  a  standard  of  210  K  rather  than  290  K,  when 
referenced  to  the  input  of  the  mixer.   In  these  units  the  system  noise 
becomes 

Tsys  =  1151   §5   =  833°K-  (B3.D 

In  addition  to  absorbing  radiation,  the  components  before  the  mixer 
emit  radiation  with  emission  coefficient  equal  to  their  absorption 
coefficient.  A  l.k   db  insertion  loss  and  29O  K  device  temperatures 
gives  rise  to  8l  K  worth  of  emission  which  must  be  subtracted  from 
the  system  noise  to  isolate  the  mixer  contribution: 

Tmix=  833  "  8l  =   753°K'  (B3'2) 

Multiplication  by  two  gives  the  standard  form  for  a  single  sideband, 

T  .   =  1506°K 

(B3.3) 
=  7.92  db. 

This  result,  when  compared  with  the  manufacturer's  quoted  noise  figure 

for  the  mixer,  indicated  that  there  is  about  a  half  db  worth  of  extra 
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insertion  loss  or  receiver  noise  in  the  system. 

To  find  the  sky  noise,  a  comparison  can  be  made  between  the 
receiver  noise  figure  and  the  effective  noise  figure  derived  from 
Equation  h]4, using  the  actual  data  from  the  experiment.  The  system 
temperature  from  the  entire  years  data  is  about  'll*00  K,  which  is  about 
250  degrees  higher  than  the  1151  degree  system  temperature  observed 
with  the  apparatus  terminated  in  a  300  degree  load.  This  additional 
noise  is  due  to  the  sum  of  the  environmental  effects  operating  on  the 
polarimeter  and  is  probably  "sky  noise"  from  weather  and  other 
atmospheric  fluctuations. 

Bh.      THE  CIRCULATOR  AND  ISOLATOR 

The  factory  specifications  of  the  circulator  and  isolator  are: 
THE  CIRCULATOR 

Ferrotec  Model  T-384L 

Serial  Number  2 

Isolation  30.3  db 

Insertion  Loss  .2  db 

> 

VSWR  1.06 

Frequency  8.5-10.1  GHz. 
THE  ISOLATOR 

E  &  M  Laboratories  Model  X112  LIA 

Serial  Number  171 

Isolation  ^lO  db 

Insertion  Loss  .8  db 

VSWR  1.05 

Frequency  9.2-9.55  GHz. 
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B5.   THE  MAGNETIC  SHIELDING 

The  magnetic  shielding  was  constructed  from  .050"  thick  Moly- 
Permalloy,  high  permeability  shielding  material  by  Allegheny  LudluJoi 
Steel  Corporation.   It  consisted  of  two  11  inch  diameter  cylinders  each 
20  inches  long  and  so  constructed  that  one  would  slide  within  the  other. 
Each  end  was  covered  with  a  spun  cap  of  the  same  material,  having  holes 
for  either  the  antenna  and  sky  horns  or  the  electronics  cables.  When 
assembled  for  use,  the  shield  formed  a  covered  cylinder  11  inches  in 
diameter  and  30  inches  long.  Moly-Permalloy,  when  properly  annealed 
has  the  following  magnetic  properties : 

Operating  Permeability  85,000-135,000 

Maximum  Permeability  300,000-600,000 

Induction  at  Maximum  Perm.        3000  Gauss 
Saturation  Induction  8000  Gauss 

Specific  Gravity  8.7^ 
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